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Abstract 

We present a QCD analysis of heavy quark mesons focussing on the B — > D* formfactor 
at zero recoil, J r o*(l). An advanced treatment of the perturbative corrections in the 
Wilsonian approach is presented. We estimate the higher-order power corrections 
to the OPE sum rule and describe a refined analysis of the nonresonant continuum 
contribution. In the framework of a model-independent approach, we show that the 
inelastic contribution in the phenomenological part of the OPE is related to the m,Q- 
dependence of the hyperfine splitting and conclude that the former is large, lowering 
the prediction for Td* (1) down to about 0.86. This likewise implies an enhanced 
yield of radial and D-wave charm excitations in semileptonic B decays and alleviates 
the problem with the inclusive yield of the wide excited states. We also apply the 
approach to the expectation values of dimension 7 and 8 local operators and to a few 
other issues in the heavy quark expansion. 
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1 Introduction 



A precision determination of the parameters of the Standard Model (SM) is mandatory for 
a stringent test of the theory and eventually to identify the effects of new physics. Decades 
of large-scale experimental efforts have been devoted to testing the SM. At high energies 
LEP has examined the gauge structure of the SM with great precision without any serious 
hint of effects beyond the SM. The flavor structure of the SM has been verified in many 
ways, and in particular at the B factories of SLAC and KEK. In both respects, the flavor 
and the gauge structure, the SM has so far passed all the tests. 

A sound theoretical foundation is needed to perform such tests along with precise data. 
Methods have been developed in heavy flavor physics over the last two decades that allow 
systematic and controllable calculations. As an example, the determination of the CKM 
matrix element V c b employing the heavy mass expansion has reached a relative theoreti- 
cal uncertainty of about two percent; its extraction from inclusive semileptonic decays is 
presently believed to be most precise. 

The determination of V c b from the exclusive decay B — > D*£u requires the knowledge 
of the form factor at the zero-recoil point where the velocities of the initial and final states 
are equal, v ■ v' = l. In the heavy quark limit the form factor is normalized to unity owing 
to the heavy quark symmetries; the precision determination of V c b requires to control the 
deviation from the symmetry limit. 

The zero-recoil version of the heavy quark sum rules was proposed back in the 1990s 
to estimate the scale of the nonperturbative corrections, which turned out significant, even 
considering that the effects are driven by the moderate charm mass. As soon as the first 
experimental B — > D*£u data were available, these early analyses suggested a value of V c b 
in agreement with the inclusive value, but were unable to make precise predictions because 
of the poor knowledge of the important heavy quark hadronic parameters. 

A model-independent accurate calculation of the formfactor may be expected from QCD 
lattice simulations provided they measure directly its deviation from unity. The value for 
the form factor obtained in existing lattice calculations leads to smaller V c b compared 
to the values from the inclusive determination. The latter meanwhile has become quite 
mature both experimentally and theoretically. 

The significant progress in our understanding of heavy mesons brought about by all 
the precision data in the B sector calls for a reappraisal of the heavy quark expansion for 
In particular, it allows us to obtain a prediction for .F(l) with an informative and 
defendable error estimate. 

In the present paper we discuss the technical details of the analysis that led to the 
results reported in Ref . [I] . Our estimate for the central value of the formfactor at v ■ v' — 1 
is appreciably lower than lattice estimates, which in turn appear marginally compatible 
with the unitarity bound. Consequently, the value of V c b extracted using our formfactor 
is larger and happens to be close to the inclusive one, although it still suffers from larger 
theory uncertainty. 

In the course of the analysis of the zero-recoil sum rule we have found a fruitful link 
among three apparently unrelated topics in heavy meson phenomenology: the value of ^(1), 
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the hyperfine splitting in B and D, and the '| > |' puzzle. In all three cases the experimental 
data are naturally described if the transition amplitudes into the 'radial' and/or D-w&ve 
states are more significant than in the naive quark models. 

The paper is divided in two parts. The first part describes the principal ingredients 
required for the evaluation of .F(l). Our general approach, outlined in Sect. |2j is essentially 
an implementation of the Heavy Quark sum rules laid down in Ref. [2J. Technically, we 
formulate it in a somewhat different way which has several advantages. To eventually 
sharpen the estimate, a number of points had to be improved. 

First, the perturbative corrections should be calculated with the hard cutoff consistent 
with the Wilsonian 'kinetic' renormalization scheme. This is a nontrivial point which caused 
confusion in the past literature. We present a method to obtain the leading a s -corrections 
together with all BLM-improvement terms as a function of the cutoff scale fi, avoiding 
an expansion in [ijvrtQ. This is the subject of Sect. [3] where the numeric results are also 
presented, including known second-order non-BLM terms. 

Section [4] discusses power corrections in the short- distance expansion of the scattering 
amplitude off the heavy quark; these correspond to the power corrections in the sum rules. 
We find a noticeable impact of the \ jw? c corrections lowering J-"(l), and estimate higher- 
order effects in the framework of the recent analysis [3] of I/tuq and 1/rriQ corrections in 
the inclusive decays. 

Section [5] presents the novel evaluation of the inelastic contributions to the unitarity 
relation, which uses as an input the hyperfine splitting of B and D mesons. We found that 
the usual nonlocal D = 3 correlators for heavy quark mesons are rather large, and this yields 
a large inelastic contribution which significantly exceeds the naive estimates of the past. 
This in turn lowers the expected central value of .F(l) down to about 0.86. We have also 
refined the estimate of the D^*'tt continuum contribution to the inelastic contribution, which 
supports its overall numerical significance. Independently of the hyperfine constraints, the 
continuum contributions lowers the unitarity upper bound for .F(l). 

The conclusions of the first part focussing on are summarized in Sect. [8j with the 
numeric values given in Eqs. ( pi] ), ( p2j ) and ( p3| ). 

Section [6] opens the second, more theoretical part, where we scrutinize, in a model- 
independent way, the higher excited heavy quark states contributing to the inelastic tran- 
sitions. In the heavy quark limit these belong to the radial or D-wave states. This is the 
framework which allowed us to relate the mass dependence of the hyperfine splitting to J-"(l) 
in an informative way. Within the same framework we could link the observed enhance- 
ment of the nonlocal effects to the significant inclusive yield of higher-excited states beyond 
D, D* and their P-wave excitations, thereby substantiating a possible resolution of the 
so-called £ | > §' puzzle. As a byproduct of the model-independent description we analyze 
the nonfactorizable effects in the higher-dimension expectation values, and we give numeric 
estimates for some representative combinations. The non-resonant D^tc continuum has 
been studied in the same heavy-quark limit setting, relying on the soft-pion approximation. 

In Section [7] we discuss certain theoretical aspects related to both the present and the 
lattice analyses of ^(1). We also briefly mention what can be gained applying our approach 
to the vector B^tD transitions. The main conclusions are summarized in Sect. [U 
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Figure 1: The analytic structure of T zr (e) and the integration contour yielding the sum 
rule. Distant cuts are shown along with the physical cut. The radius of the circle is Em- 

The Appendices contain several of the details omitted from the main text. They mostly 
concern the perturbative calculations; there we discuss some of the conceptual aspects and 
also give the concrete expressions required in the analysis of ^(1). 

2 The framework 

We consider the zero-recoil (q = 0) forward scattering amplitude T zr (e) of the flavor-changing 
axial current C7756 off a B meson at rest: 

T z \s) = fd 3 x J dx e -iMM B -M D *- £ )_l_ {Bl i iT c m5 b(x) hmc(0)\B) , (1) 

where e is the excitation energy above Mo* in the B — > X c transition (the point e = 
corresponds to the elastic B — > D* transition). The amplitude T zr (e) is an analytic function 
of e and has a physical decay cut at e> 0, and other 'distant' singularities at |e| >2m c . The 
analytic structure of T zr (e) is illustrated by Fig. [TJ 
We consider the contour integral 



Io(e M ) = -— <t T™(e)de (2) 




with the contour running counterclockwise from the upper side of the cut, see Fig. [T] Using 
the analytic properties of T zr (e) the integration contour can be shrunk onto the decay cut; 
the discontinuity there is related to the weak transition amplitude squared of the axial 
current into the final charm state with mass Mx = Mp* +e. If we explicitly single out the 
elastic transition contribution B^-D* at e = then 

1 f £M 

I (e M )=F 2 (l) + l mel (e M )i Imei(e M ) = ^— / discT zr (e)d £ (3) 

holds, where iinei(£itf) is related to the sum of the differential decay probabilities, in the 
zero recoil kinematics, into the excited states with mass up to Md*+Em- 
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The OPE allows us to calculate the amplitude in (UJ) - and hence Iq{sm) - in the short- 
distance expansion provided |e| is sufficiently large compared to the ordinary hadronic mass 
scale. It should be noted that strong interaction corrections are driven not only by |e|, but 
also by the proximity to distant singularities. Therefore, Em cannot be taken too large 
either, and the hierarchy e M <C 2m c has to be observed. 

The sum rule Eq. (|3| can be cast in the form 



= ^/He M )-hn^e M ) (4) 

which is the master identity for the considerations to follow. Since I\ ne \{sM) is strictly 
positive, we get an upper bound on the formfactor 

•F(l) < Vh(e M ) (5) 

which relies only on the OPE calculation of Iq. Note that this bound depends on the 
parameter em, while Eq. Q is independent of Em since the dependence in I and Ji ne i 
cancel. Furthermore, including an estimate of I- me i(^M) we obtain an evaluation of J-"(l). 

The correlator in Eq. ([T| can be computed using the OPE, resulting in an expansion of 
T zr (e) in inverse powers of the masses m c and m^. This yields the corresponding expansion 
of Iq(em)- This OPE takes the following general form: 

t V"^ 2M {■B\Ok\B)[ 1 

Io(e M ) = £T (£M,fj) + V C k (e M ,fj) — — j- 3 (6) 

k m Q 

— i^ n { £ M,y) - A.!/ m 2 {em. A*) - A!/ m 3 (e M , A*) - A x / m 4 /x) - ... 
= e A ert (e M ,f,)-A A (£ M ,v), 

where Ok are local 6-quark operators b...b of increasing dimension dk > 5, Ck{EM,n) are 
Wilson coefficients for power-suppressed terms, and £^ er is the short-distance renormaliza- 
tion (corresponding to the Wilson coefficient of the unit operator), which is unity at tree 
level. We have also introduced a Wilsonian cutoff fi used to separate long and short dis- 
tances. The complete result does not depend on /i since the /i-dependence cancels between 
the Wilson coefficients and the matrix elements of the operators. At tree level A^ does not 
depend on em- The choice of /i is subject to the same general constraints as that of Em, 
and therefore we will often set /i — Em, in which case we will also use ^ ert (A*) = ^A^it 1 , A*)- 



3 Perturbative corrections 

The leading perturbative renormalization factor ^ ert (e M ,/i) can be expanded in powers of 
a s . In the Wilsonian OPE all infrared physics is removed from perturbative corrections; 
the perturbative series for £^ ert is then free from infrared renormalons. The exact form 
of the perturbative coefficients depends on the concrete definition of the higher-dimension 
operators used in the OPE; we consistently assume the scheme of Refs. jH [5] often referred 
to as "kinetic" (or Small Velocity); see also Ref. j6]. Here we describe a compact method 
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to calculate £^ ert in this scheme; the details as well as the justification of the method can 
be found in Appendix |B] 

For the theoretical analysis we need to keep \x and Em distinct. Their role is different: 
Em specifies the observable Iq{em) under study, the energy integral of the spectral density. 
Iq(em) depends on e m , and the same applies to the Wilson coefficient £^ ert for its leading 
term in the OPE. On the contrary, \i is a technical tool employed in the OPE to separate 
'soft' and 'hard' physics. Since Iq(em) should not depend on /i, the same applies to the 
right-hand side of Eq. (JgJ) . While £^ ert does depend on /x once the perturbative corrections 
are included, this dependence is canceled by the /i-dependence of the power-suppressed 
expectation values. 

The ^/-dependence of J is governed by the inelastic spectral density Wi ne \: 

dIo( £ M) dJi ne i(e^) 

— H = — J~ c = W inel(£Mj- (7) 

The same holds for the ^-dependence of (em, aO as long as only the perturbative part 
of the inelastic spectral density is considered]^] 

der(£M ^W ert (e M ), (8) 



ds 



M 



where w peit (E) = w?* el (e) 



Although both em and fx should be sufficiently large, in a fixed-order perturbative 
expression one can formally set them equal to zero in ^ ert (£jn, fi). It is evident from the 
definition of Iq{em) that this would yield a 'purely perturbative' renormalization factor t]a 
for the zero-recoil axial current to this order, as it is commonly defined in HQET: 



er^o) = v i ( 9 ) 



We need, however to keep both mass parameters nonvanishing. Eq. (|8J) allows one to pass 
to non- vanishing Em- In general, the /^-independence of the overall result is used to fix £^ ert 
at non-vanishing \i. 

A relation similar to Eq. ([£) holds for a physical choice of the arguments in ^ ert if 
the perturbative calculation of the renormalization factor tja is performed with an infrared 
cutoff /i, and provided /i is close to Em- There is also a perturbative contribution to 7 ine i, 
whose soft part is removed by the same cutoff. The precise relation is not trivial, however, 
even at one loop. An extensive discussion is given in Appendix [B] 

A direct calculation of the Wilson coefficients implies an infrared cutoff on the internal 
gluon momentum in the Feynman diagrams. The contribution of soft gluons is subtracted 
from 7] a by a term Srj 8 ^ ft . In one-loop diagrams, see Fig. [2] the kinetic scheme cuts off the 
Feynman integral at \k \ </i yet the integration over ko is always performed from —00 to 00. 
While the integral is dominated by small ko<\k\ , a power-suppressed contribution comes 



1 This is true in general for e M >• [/,, or if we neglect 0{a s ) contributions to the Wilson coefficients of 
the power suppressed operators. 
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Figure 2: a) Vertex diagram contributing to t]a- b) poles in the complex ko-plane encountered in 
the one-loop diagrams. The standard Feynman propagators lead to contour a (blue); upon closing 
it into the lower half-plane the integral equals the sum of the three residues. The calculation of 
St]^a (/i) requires to pick up only the 'near' gluon pole (red), and is obtained by integration over 
contour j3 (the bypass in green). 

also from large fco~ m Q- This piece is not soft and would not be properly accounted for if 
attributed to the matrix elements in the OPE. 

This brings us to an important point: the subtraction of the full \k\ < fi one-loop 
contribution to tja would not yield the correct Wilson coefficient, leading to spurious terms, 
formally of order 0([i 3 /tuq). The correct subtraction includes only the residue of the 
'near' gluon pole at ko = \k\—iO in the calculation of the integral over k Q . The other two 
residues encountered in the conventional Feynman diagram calculation upon closing the fco 
contour into the lower half-plane, see Fig. [2}d, correspond in fact to hard physics with gluon 
virtuality ~2mg and must be left to the hard Wilson coefficients. 

The different physics associated with the 'distant' poles in the diagram located near 
k = 2rriQ + k 2 /2rriQ can be intuitively understood: picking up the corresponding residue 
would leave the gluon propagator hard, k 2 ~ 4mg for small k. These poles are actually 
related to the divergence of the power expansion in k/rriQ when the soft scale \i is increased 
towards tuq. They contribute terms ~ d 3 fc/mg oc /t 3 /mg and are seen starting at order 
1/ttiq in the OPE. A detailed discussion of this technical point is given in Appendix p| The 
validity of this prescription in the kinetic scheme is explicitly checked by the /^-independence 
of the OPE relations beyond order 1/wq. The matching of the /^-dependence of £4 to that 
of the OPE expectation values is demonstrated in Appendix [Cj We have also verified this 
for the 1 / rriQ expansion of the heavy hadron mass. 

As a result, the Feynman integration must be modified to exclude the residues of the 
distant poles in the ko plane from the calculation of 8r^ . The complete one- loop expression 
for ^ ert (eM)rf accounting for the e^-dependence, Eq. (J8J), takes the following form: 

e A eit (e M , //) = V l - 2 5 V f\ri - [" de w^(e) } (10) 
where // corresponds to the excitation energy of the inelastic transition with emission of a 
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X X 



Figure 3: The poles in the integrand and the integration contour in the complex ko plane for 
calculating the Wilson coefficient. At k 2 + \ 2 > li 2 the integration is performed in the standard 
way, dashed line. At k 2j t-X 2 <Li 2 the 'near' gluon pole (red) is moved above the real axis, and the 
integration contour changes to pass below it, along the green short-dashed path. A stands for the 
gluon mass. 



gluon with energy LU — \k\—fJ,, 



/ i ' = / i + V / ' m c+/ i2_m c- (11) 



The last term in Eq. (10) describes, at fi — EM, the recoil correction in the relation between 



the normalization point and the upper limit of energy integration in Iq. Let us note that 



the recipe for calculating the combined contribution of the first two terms in Eq. ( 10 ) can 
also be formulated as changing, for \k\< ft, the bypass prescription for the gluon pole at 
k = \k\ from k +i0 to k — zO, see Fig. |3| In other words, the proper Wilsonian prescription 
consists in discarding the gluon pole contribution for gluons softer than the cutoff. 
The one-loop renormalization factor t]a without a cutoff is well known: 



VA = 1+ \CA ( In ^ - r , + (a 2 ) ; (12) 



4 7r \mj- m c m, 
the explicit calculations yield for 5r] s ^ ft (fi) 

5r£*tu)= Cf9 * f d3 ^ ( 1 | 2 | 1 V gjgwgfi | 2 | 1 

1 " 4 (27r) 3 |fc| V m c 3m c m b m 2 ) vr 4 \m 2 3m c m b mf 

(13) 

and the 0{pc s ) inelastic spectral density is 

W ^(e) = C F ^ \ (2M 2 + 3m 2 b + 2m b m c + m 2 c ) , M = m c +e; (14) 

ir 12 M 6 mi 

M is the invariant hadronic mass in the final state. The integral of the one-loop spectral 



function that appears in (10) is 



= { ^%^±^ [2 £M ( £M+ 2m c ) (15) 

J , ir I 2Am z b (e M +m c ) 2 L 



mz 



+ 



(ml — 3ml — 2m b m c + 4m cj u' + 2fs! ) -| (3m b — m c )(m b +m c ) ^ m c +fi' 

12 m? m c +e M 



(m c + fx') 
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Combining these in Eq. (10) we obtain the one-loop £^ ert (e m , y) as an explicit function of 
n/rriQ. The numeric dependence on [i can be seen in the first-order plot in Fig. |4j 

Our analysis of the one-loop corrections can be readily extended to include higher- 
order BLM corrections which describe the effect of the running of a s in one-loop diagrams; 
a complete BLM-summation is also possible. A detailed discussion of the technique, in 
the context of the Wilsonian OPE, can be found in the Appendix of Ref. [7j{^] We have 
recapitulated the salient points in Appendix |Aj In practice, a fictitious gluon mass A is 
introduced in the one-loop diagrams, and eventually a weighted integral over A 2 is taken. 
Consequently, the above discussion concerning the Wilsonian cutoff and the OPE applies 
also to the BLM corrections of any order. 

The case of a massive gluon requires the following kinematic modifications. Since the 
gluon energy is now 



uj = \Jk 2 + X 2 , 

the range of soft gluon momenta k shrinks, \k\ < ji 2 — X 2 and the cutoff in the diagrams 
is now triggered by 

9(fi 2 -X 2 -k 2 ). (16) 
Notably, no subtraction is necessary at A > fi. The recoil energy is also modified and \J in 



Eqs. (11) is given by 

// = \Jm 2 +ii 2 — A 2 — m c +fi; (17) 



the perturbative inelastic spectral density is given in Appendix [C) Eq. ( |A~13)). 



In the Feynman diagrams the relevant gluon pole is now located at ko = y k 2 + X 2 rather 
than at k = \k\; this makes the explicit expression for 5f]^ more cumbersome: 



8r,T^) = C F g 2 s f 

Jkn 



d 3 k f2m 2 -2m b k -2k 2 + X' 2 



k0<fl (2^2k V (2m b k -X 2 ) 2 

2m 2 - 2m c k - 2 k 2 + A 2 2m b m c - (m b +m c ) k + § k 2 + A 2 \ , . 
+ (2m c k -\ 2 ) 2 (2m b k -X 2 )(2m c k -X 2 ) V ^ } 

the integral can be solved analytically, resulting in a lengthy expression. The expression 



for rjA at nonvanishing gluon mass is given in Eq. (A. 14) of Appendix 



Combining these elements in Eq. (10) we obtain the one-loop correction to £^ crt at 
arbitrary gluon mass A, [sm-, A<; A 2 ). Assuming /j, — e M , for A 2 >/z 2 the last two terms in 
Eq. (10) are absent at A 2 >/i 2 and one simply has ^ e (//; A 2 )=?7^(A 2 ). Using the formulas 



of Appendix [A] one readily obtains the BLM corrections of arbitrary order or the resummed 
result. We performed the final integration over the gluon mass numerically. 

The explicit expression for ^(A 2 ) at small A 2 shows non-analytic terms in A 2 , starting 
with \- In A 2 . They signal the infrared sensitivity of 714 and the emergence of infrared 

m Q 



2 There was a typo in Eq. (A. 6) of paper [8] for the BLM-resummed expression which unfortunately 
propagated into the later paper [7], Eq. (A. 20); in that equation Ay must be replaced by Aq CD (the 
conventional MS one) in the denominator of the power term. The correct expression is given here in 
Appendix [Aj Eq. ( |a3| ). 
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renormalons, which in turn makes it impossible to assign a definite value to the purely 
perturbative va, and leads to a significant numerical instability from the higher-order cor- 
rections. Non-analytic terms are also present in Srj^ ft , and they precisely offset those in 



^(A 2 ): the combined £^ er (e^f,//) in Eq. (10) is an analytic function of A 2 in the vicinity of 
zero at any positive \x. The radius of convergence of the Taylor series in A 2 is precisely /i 2 . 
The cancellation of all the non-analytic pieces is most simply seen directly in the integral 
representation for tja — ^^(/x), upon closing the contour over the residues at positive k . 
The recoil integral describes the effect of shifting the normalization point and is purely 
short-distance. Indeed, the one-loop spectral density w pcrt (e; A 2 ) is explicitly analytic at 
small A 2 for e> fi where \k\ is of order /i. Consequently, there are no formal obstacles in 
deriving either any higher-order BLM coefficient or the fully resummed BLM value for the 
Wilson coefficient ^ ert (^Af , A*); apart from the ultraviolet domain, its perturbative series 
has a finite radius of convergence, at least in the BLM approximation. 

The advantage of the method described in this section is that it allows to calculate the 
full /i-dependence of the Wilson coefficient. It does not apply to non-BLM corrections, 
starting with 0(a 2 s ). In that case, the /i-dependence of £^ ert has to be determined, order 
by order in 1/mg, applying the normalization point independence of the OPE relations, 
and using the initial condition Eq. (J9|). This was done for the non-BLM 0(a 2 ) corrections 
in Ref. [9], through 0(l/rag); in particular, the two-loop spectral density w pert (e) was 
calculated there to this accuracy. The corresponding non-BLM coefficient was found to be 
small numerically, which suggests that omitted terms 0(a 2 /i 3 /mg) and higher should not 
have a noticeable impact. 

3.1 Numerical analysis 

The perturbative corrections i A ert {sM) appear to be small for values of Em between 0.6 GeV 
and 1 GeV. Taking, for instance, Em — A 4 — 0.75 GeV, m c = 1.2 GeV, nib = 4.6 GeV and 
assuming a s {rrib) = 0.22 we get 



£P° rt = 1 - 0.022 + (0.005 - 0.004) + 0.002 - 0.0015 + ... (19) 

Here the first term is the tree-level value, the second is the 0(a s ) term evaluated with 
a s = 0.3, which corresponds to the strong coupling evaluated at an intermediate scale 
between m c and m&. The two values in brackets show the shift, relative the one-loop 
evaluation with a s = 0.3, due to 0(a 2 s ) corrections (positive for the BLM part and negative 
from the non-BLM contribution); the last two terms are the 0(/5 2 af) and 0((3qQz^), which 
may serve as an estimate of even higher-order perturbative corrections. 

Fig. j^J shows the dependence of \J £^ crt (/i) on fi at different orders in a s assuming 

a^ s (m fe ) = 0.22. For \i between 0.7 and 0.8 GeV the value of J ^ ert (/i) is close to 0.98, and 
we associate to it a rather conservative 1% uncertainty: 



^ ert (0.75GeV) = 0.98 ±0.01 
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Figure 4: J^ifi) as a function of fi for m b = 4.6 GeV, m c = 1.2 GeV, a s (m b ) = 0.22. The 
curves represent the one-loop result evaluated with a s = 0.3 (blue), one-loop plus first order BLM 
(green), complete O(a^) (red), two- loop plus third-order BLM (maroon). 

We emphasize that the observed stability of the perturbative expansion applies only to the 
perturbative renormalization factor in the Wilsonian OPE, and that the quoted value refers 
to the specific renormalization scheme which is used in the present analysis. 

Let us note here an observation that is relevant when the terms 0{l/m%) and higher are 
accounted for. As mentioned in the previous subsection, the calculation of the perturbative 
factor £^ er (/x) which neglects the subtleties related to the 'distant' poles yields the correct 
result through order /i 2 /mg, and one needs an improved method only to account for the 
terms /i 3 /mg and higher. Numerically the two methods yield very close results to order a s 
and /3o«s, but the difference starts to accumulate systematically from order (3$^ on, even 
though parametrically the difference between the two methods is still only of order l/m%. 
For instance, the fully resummed results for the two prescriptions to calculate £^ ert (/x) differ 
by about 2%. This is in contrast with the case of the total b^rciv width, where there is 
no visible difference between the result truncated at order a 2 s and the fully resummed one. 
The reason is that the current renormalization at zero recoil has an intrinsically lower scale 
driven by the charm mass, in contrast to nib for r s \(B). 




4 Power corrections to Iq 

In this section we investigate the power corrections in the right-hand (OPE) side of the 
sum rule Eq. The leading power corrections to I were calculated in Refs. [10J |2] to 
order I/tuq and in Ref. [TT] to order 1/rriQ and read 

+ / j_ + 2 + ^ 

3m 2 4 \ m c 3m c mb m h 



,3 1 „3 



Pd - iPls 1/1 1 3 . 3 



Al/m3 = Ami + 12^ W + ^ + (PD + PLS) ' (21) 
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The nonperturbative parameters [Mq, p|> and p\ s all depend on the hard Wilsonian 
cutoff. In the renormalization scheme we have adopted the inequalities fi^ip) > /^(/-t), 
Pd(aO > — Pis(p) hold at arbitrary normalization point /i. The nonperturbative contribu- 



tions in Eqs. (20), (21) are therefore positive. 

To remain on the conservative side for numeric estimates we can adopt the low values 
j4(0. 75 GeV) = 0.4 GeV 2 , p^(0.75 GeV) = 0.15 GeV 3 , and the quark masses m c = 1.2 GeV 
and rrib = 4.6 GeV (the scale dependence of the latter plays a role here only at the level 
formally beyond the accuracy of the calculation). The dependence on p? G and on p\ s is 
minimal and their precise values do not matter; we use for them 0.3 GeV 2 and -0.12 GeV 3 , 
respectively. We then get 

A 1/m2 = 0.091, A 1/m 3 = 0.028. (22) 

It is interesting to compare these estimates with those derived from the constraints 
on the expectation values of dimension 5 and 6 operators from the semileptonic S-decay 
moments. If we employ the values of the OPE parameters extracted from the latest of- 
ficial HFAG fit to inclusive semileptonic and radiative decay distributions [121 for 
p = 0.75 GeV we find 

A 1/m 2 + A 1/m s = 0.102 ± 0.017, 



which is consistent with (22). As discussed in [T3], this HFAG fit to semileptonic moments 
depends on several assumptions and does not yet incorporate certain higher-order effects 
that may be important, including the complete a 2 -corrections [H]. In particular, more 
realistic ansatze for the theoretical correlations have been considered in [15J, leading to 
larger values of the A A , with bigger errors. Typically, one then has 

A 1/m 2 + A 1/m3 = 0.11 ±0.03. (23) 

On the other hand, combining the semileptonic moments alone with a high precision deter- 
mination of the charm mass [H] yields 



A 1/m2 = 0.090 ±0.013, A 1/m3 = 0.029 ±0.008, (24) 



in remarkable agreement with Eq. (22). 



The important question is how well the power expansion for the sum rule converges. 
Recently, the OPE for the semileptonic S-meson structure functions has been extended to 
order 1 / nig and 1 / tjiq [TTJ [3] . Applying the analysis to the structure functions mediated 



12 



by the axial and by the vector currents separately, we fincj^] 
16m 4 A 1/m4 = _(3+§y+2/ 4 )mi + (ly~ly 2 ~ly 3 -2y 4 )(m 2 +m 5 ) - (1 + %y + |y 4 )m 4 

+ (l + |2/-2/ 4 )^6 - |?/^7 + (l + |y-|2/ 4 )^8, 

16m* A 1/m a = -(2+|t/-| 2/ 2 )r 1 + (6+% + ^ 2 + §y 3 + 4 y 4 +2y 5 )(r 2 -r 3 ) 
+(6+4y-f y 2 + | y 3 + 2y 5 )r4 - (2+% + | 2 / 3 + | 2/ 4 +y 5 )r 5 
-(2+fy+f y 2 + |y 3 + 2/ 4 +y 5 )r 6 + (2 + ly + ly 2 + ly 3 +y 4 + y 5 )r 7 
+(|-| 2 /+|y 2 )r 8 - (2+y + ^ 2 -fy 3 -§y 4 -3y 5 )(r 9 -r 12 ) 
-(2+y+^ 2 -|y 3 -| 2/ 4 -3 2 / 5 )(r 10 -r 11 )-(2+6y+fy 2 -f 2 / 3 -|y 4 -% 5 )r 13 
-(2-% + fy 2 -| 2/ 3 -2 2/ 5 )r 14 + (|-| 2/ +% 2 -| 2/ 3 -| 2/ 4 -3 2/ 5 )r 15 
+(§ + ^+V-§y 3 -fy 4 -2 2/ 5 )r 16 + (l-|2/-f2/ 2 -|2/ 3 -2/ 4 -2/ 5 )n 7 
-(| + l2/+¥2/ 2 -|2/ 3 -|2/ 4 -32/ 5 )^ 8? ( 2 5) 

where y = m c /mb and the D = 7 and D = 8 expectation values mi_g and ri_i§ are defined 
in Ref. [SJ. These can be evaluated in the ground-state factorization approximation. Using 
the expressions given in Ref. [3J we obtain the estimates 

A 1/m 4 ~ -0.023 , A 1/m5 ~ -0.013. (26) 

It is worth noting that retaining only the terms suppressed by the powers of l/m c (i.e., 
evaluating the higher-order corrections in the limit m& — > oo) yields a perfect numeric 
approximation to the full expression. In the ground state saturation approximation the 
dominant contributions to A]/ m 4 and A]/ m 5 are those of m 48 and r 2) io, respectively, without 
significant cancellations. We then observe that the power series for J appears well-behaved 
at the required level of precision. 

For what concerns the loop corrections to A A , the a s -correction to the Wilson coefficient 



for the kinetic operator in Eq. ( 20 ) was calculated in Ref. [H] and turned out numerically in- 
significant. Generally larger a^-corrections are expected in the chromomagnetic and Darwin 
coefficient functions. However, the dependence on \Xq in Eq. (20) turns out negligible; there- 



fore perturbative corrections are not expected to introduce significant numerical changes in 
the estimate of Ai/ m 2. At order l/m%, even if radiative corrections change the coefficient 
for the Darwin term by 30% the effect on the sum rule would still be small. 

Taking into account all the available information, our estimate for the total power 
correction at £M = 0.75GeV is 

A A = 0.105 (27) 

with a 0.015 uncertainty due to higher orders. On the theoretical grounds, larger values of 
/^ 2 and/or p 3 D are actually favored; they tend to increase A A . Combining the above with 
the perturbative corrections we arrive at an estimate for J and, according to Eq. ([5]), at a 
bound on the form factor, which in terms of the central values at = 0.75 GeV is 

^(1) < 0.925. (28) 



3 We thank S. Turczyk for providing us with the input needed for this calculation. 
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As stated above, the upper bound in Eq. ^ depends on em, becoming stronger for 
smaller em- It is advantageous to choose the minimal value of Em for which the OPE-based 
short-distance expansion of the integral (|2| for Iq(em) sets in. This directly depends on 
how low one can push the renormalization scale fi while still observing the expectation 
values actual /^-dependence in the kinetic scheme approximated by the perturbative one. 
Since in this scheme /^-(/f) > p^{p) holds fo r arbitrary /i, in essence this boils down to 
the question at which scale /i m i n the spin sum rule and the one for \Xq get approximately 
saturated, e.g. /i G (/i m in) — 0.3GeV 2 . The only vital assumption in the analysis is that the 
onset of the short-distance regime is not unexpectedly delayed in actual QCD and hence 
does not require Em > 1 GeV. This principal question can and should be verified on the 
lattice. This will complement already available evidence from preliminary lattice data [TS] 
as well as from the successful experimental confirmation [19] in nonleptonic B decays of 
the predicted | -dominance. 



5 Estimates of ii ne i 

We now turn to the actual estimate for the inelastic contribution. On general grounds [10J 
ii ne i is expected to be comparable to the power correction A A considered above; therefore 
the inelastic contributions should be important numerically. 

Our starting point is the first moment of the scattering amplitude spectral density given 
by the contour integral 

1 



h(l*) = -7r-. f T^(e)ede; (29) 
2m J 

\e\ = £M 

we can write 

Wew) = , (30) 

where e is an average excitation energy which depends on Em- For moderate Em the integral 
is expected to be dominated by the lowest 'radial' excitation^] of the ground state, with 
£ ~ £md ~ 700 MeV. The first moment Ii(em) can also be calculated in the OPE [2]; the 
result reads 

Il _ -(pIg + Pa) , '^PL-Pip , pL+pIg+Ps+Pa ( 1 , 2 | 1 \ ( 1 \ 
?>m? c 2>m c mb 4 \ m c Sm c mb my \ m Q J 

The nonlocal zero momentum transfer correlators p^, p^ G , p| and p\ have been introduced 
in [2] and are given by 

PL = J d 4 x ^-(B\iT&if*b(x), b^bm\B)\ 



4 These excited states play an important role in the analysis of power corrections in the HQE, and we 
clarify our terminology. In the heavy quark limit these are either the true radial excitations of the ground 
state, or the counterparts of the I?- waves. The former have spin-parity of the light cloud \ + while for the 

latter it is | + . At finite quark masses by "radial excitation" we refer to the descendants of any hyperfine 
multiplet member of these heavy-quark states. More details are addressed in Sect.[6j 
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^PsMkJ + IpaWjI - S u 6 jk ) = [ d 4 x -j-iBltTibaMix), 6o- J -S l 6(0)}|S)', (32) 

where 5 denotes the chrome-magnetic field strength operator. The prime indicates that 
the ground-state contribution is subtracted - otherwise the integral diverges at large xq 
where the correlators approach a constant determined by the ground-state factorization 
contribution. 



At higher orders in 1/toq the expansion (31) will include, along with the local expecta- 



tion values of higher dimensional operators, more intricate nonlocal T-products. The latter 
are poorly known. Moreover, the quantum-mechanical interpretation of these relations tells 
us that there will be significant cancellations among different terms in higher orders; say, 
the factorizable terms must drop out. 

Since our goal is only to obtain a reasonable estimate, we discard higher-order corrections 
and keep only the leading 0(1/171%) terms. This implies that the expectation values can be 
considered in the static theory. The nonlocal correlators appear in I\ because the energy 
variable e is defined with respect to the physical threshold go = Mb — Mo* rather than 
relative to the parton- level threshold at qo — mb—m c native to the OPE. The terms given 
by the local operators cancel out in I± to the leading order, as becomes transparent in the 
quantum-mechanical interpretation of Ref. [2J, Sect. 6. The latter also explains the explicit 
field-theoretic expression for w- me \(e), which is given in Sec. 6.1. It is important to stress 



that the third term in (31) is positive since the combination p^ n + p^ G + pg+ p\ is actually 
equivalent to the correlator of two identical operators b(aif) 2 b. 

An important piece of information is provided by the heavy quark mass dependence of 
the hyperfine splitting AM 2 , which allows us to estimate the overall magnitude of the non- 



local correlators in Eqs. (31). The 1/rriQ scaling of the vector-pseudoscalar mass splitting 
is nonperturbatively affected by the values of two D = 3 spin-triplet parameters, p\ s and 



pIg+Pa- 



Mb* - M B = \%L + p «° + f\ Pls + O ( -1 ) , (33) 
Suit, 6mi \ m t J 

and likewise for charm. Since the spin-orbit expectation value is reasonably constrained by 



the heavy quark sum rules, (33) yields information on the combination 

-(pIg+Pa)- 

A preliminary analysis was outlined in Ref. j6] and indicated a large value exceeding the 
naive expectations. We reconsider it carefully in Sect. [6] and confirm the observation. This 



combination enters directly the expression (31) for 7 ine i and is particularly important, as 



will become clear in the following subsection, where we combine the theoretical expressions 
for Jinei with the numerical analysis of the hyperfine splitting to arrive at our estimate for 



Jinei. Additional technical details are given in Sect. 6.3 
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5.1 Numerical estimate of Ii ne i 

It turns out advantageous to analyze 7i ne i from the perspective of the BPS approximation 
for B mesons [20]. In the BPS limit p^ — p G and, since p 2 —p G = ((an) 2 ) in the B meson, 
one concludes that b(a7r)b\B) =0. This would also imply —p\s = Pb- The deviation from 
the BPS limit is quantified by the smallness of the difference p 2 —p G compared to p 2 itself 
[20J. Many remarkable relations hold in the BPS limit; for instance, among the spectral 
densities of the correlators of M 2 b and baBb that we will introduce in Sec. 6, 

pf\u) = p [ p(u) = pf\u), pf\u) = p%\u) = pf\"\ (34) 
and among their integrals determining the p 3 correlators: 

pl G = -2pl n , pl G +p\ = -(pL+p 3 s ). (35) 



In the BPS limit I\ in Eq. (31) is then given by the same combination of the nonlocal 



correlators that drives, besides p 3 LS , the hyperfine splitting to order 1/rriQ, cf. Eq. (33): 



The second term in Eq. (31) for I\ is of the first order in the deviation from BPS; as 
such it is not sign-definite. However, it is suppressed by the 6-quark mass. The third term 
is of the second order in the BPS violation and is positive; it comes with a large coefficient. 
Therefore, the full expression develops only a shallow minimum where the whole sum differs 



from the BPS value by a factor of no less than 0.93, see Section 6.3.1 In fact, Ii may exceed 
the BPS value by a larger amount, although our analysis favors a negative sign for the second 
term. This typically results, at small deviations, in an overall slight decrease. 

A simple minimal - and most physical - ansatz for the spectral densities determining 
the correlators elucidates the role of the constraints the correlators obey to. It assumes 
that they are saturated by a single multiplet of excited states, for each angular momentum 
of light degrees of freedom. Apart from the excitation mass gap the relevant contributions 



are then determined by three residues; they are introduced in Sect. |6.1| and are denoted by 
P, G (for the radially excited | + ) and by g (for the | + state). In the BPS limit P = G. 



At a fixed hyperfine constraint on — (p^ g + Pa) the full expression for I\ in Eq. (31) 



depends on two dimensionless ratios, P/ G and on the relative contribution of the | state 
proportional to g 2 . The minimum value for 1\ occurs where the latter vanishes, g = 0; 
whenever | + dominates, I\ uniformly approaches its BPS value. The value of P/G for 



which there is a minimum depends only on the ratio of the quark masses, see Eq. (85). 



Fig. pi shows the variation of Ii/lf PS with P/G for a few values of the relative contribution 



v of the | state into the combination determining the hyperfine splitting. More details of 



the analysis are given in Sect. 6.3.1 



Neglecting a possible few percent relative decrease in I\ we should, therefore, adopt the 



BPS relation Eq. (36) as a reasonably accurate lower bound estimate; this leads to the 
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r / n ^ rr , , n ^ 0.45 GeV 3 + (/«+0.2)- 0.35 GeV 3 

W^f- 0.75 GeV) > V 3 ^ 2 1 , (37) 



where the dimensionless k, introduced in Eq. (73), parameterizes the exact value of the 
correlator — (p^+p^)- ft is uncertain due to unknown higher power corrections and due to 
the limited accuracy of the perturbative renormalization. Here Em is assumed to be around 



0.75 GeV to include the families of the lowest 'radial' excitations. Our analysis in Sect. |6.2 
suggests that k is relatively small, between —0.4 and 0. 



Equating e in Eq. (37) with e ra d — 700 MeV we estimate 

/incl>0.14. (38) 

We recall that, in contrast to the OPE for Io(fi), this estimate assumes only the leading 
Phadr/^o pieces in the transition amplitudes. The subleading l/m c corrections can be 



numerically significant - this is illustrated, for instance, by Eq. (22) - and can potentially 
modify the actual / ine i by as much as 30% of the estimate, even though the inclusive sums 
like Jinei are usually affected less than the individual transitions. 

The precise value of /i 2 is not yet well known; it mainly affects the degree of proximity 
of actual B mesons to the BPS regime. The BPS expansion would become more qualita- 
tive than quantitative if /z 2 eventually exceeds 0.45 GeV 2 by a significant amount. At the 
same time, as illustrated in this section, this would not affect significantly our estimates. 
Moreover, larger /i 2 lowers the model-independent upper bound which only assumes posi- 
tivity of the inelastic contribution. Complement arily, from the full set of the heavy quark 
sum rules we should expect larger transition amplitudes to the excited states at larger /i 2 . 
This conforms the physical intuition which suggests to quantify the mass scale Phadr 
governing the strength of the suppressed transition amplitudes oc fiha,dr/ m Q- There is no a 
priori reason to have small power corrections in .F(l) at large /z 2 . 
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5.2 Continuum D^ir contribution 



The resonant states are expected to dominate the inelastic transitions at low excitation 
energy. Continuum effects are formally 1/N C suppressed and are usually numerically small, 
unless the chiral singularity for the soft pion is strong enough; in the case of J ine i it is 
only logarithmic. Here we give the continuum states D^n a dedicated consideration since 
they populate the lowest energy domain and are characterized by an average excitation gap 
that can be noticeably lower than e ra d- Moreover, these states contribute to the deviation 
from the BPS regime, possibly dominating the deviation in the low-energy regime. In this 
subsection we compute their contribution to the zero recoil sum rule with the soft pion 
technique (see [21] for a review). Complementary theoretical considerations can be found 



in Sect. 6.4 Here we follow and extend the analysis of Ref. [TO] . 

Both Dtt and D*tt channels contribute to 7i ne i. The Dir amplitude is given by the sum of 
the two pole graphs in Fig. [6} The pion vertex is parameterized by the effective Lagrangian 

C x = 2g D , Dn (M D D;Dd"7T +rM B B;Bd fM 7r). (39) 

where r = gB*B-n/ 9d*d-k'i heavy quark symmetry implies r = l. In the heavy quark limit the 
two diagrams in Fig. [6] cancel each other; all inelastic transitions vanish due to heavy quark 
symmetry. A nonvanishing result emerges once the mass shifts in the virtual propagators 
of heavy mesons are accounted for, or due to r^l. The amplitude, for a charged pion, 
then becomes 

1 (D-n+\A\B+) = -g D « Dn pJ- T —) , (40) 



where (E n ,p n ) is the pion four-momentum and we have neglected subleading terms. We 
have set the weak vertex to unity, which is legitimate to the first order in the devia- 
tions of .F(l) from unity. In terms of the D meson energy, Ed, the excitation energy 
is e = Ed + E n — Md* — E^ + Mq — Md* + p%/2Mr>. A represents a power correction: 
A ~ 2/3}j,q (l/m c +l/rnb) +O(p%/mo). Expanding the amplitude in powers of I/tuq we 
would get, for r = l, 

' (D-7i + \A\B + ) = -g D « Dn p 7T A 



2 ' 



2yJM B M D e 

which has the correct scaling oc I/tuq for non-diagonal transitions that violate the heavy 
quark symmetry. However, the hadronic mass scale in the denominator of the amplitude 
is peculiar: it is the pion energy E n , which for a light pion can be significantly lower than 
the typical QCD mass gap e ra d. 



The amplitude of Eq. (40) at r — 1 gives rise to the spectral density 



mel 12tt 2 lPwl e 2 (e + A)2' l4iJ 

the logarithmic chiral singularity that emerges upon integration in the heavy quark limit 
is regulated by the spin symmetry breaking term 2/3 (l/m c +l/m.&) in A, even for a 
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Figure 6: Pole diagrams for the B^-D+ir amplitude. 



massless pion. As noted in [TO], however, the constant term is numerically larger than the 
chiral logarithm, In A. The contribution due to neutral pions is half of the one related to 
charged pions. 



The B^rDit amplitude in Eq. (40) in fact receives additional relativistic corrections; it 
is unique because of the soft pion enhancement we have discussed. Since this enhancement 
is mild in the integrated probability, the regular contribution to the amplitude should be 
included. In the actual calculations we use the complete relativistic propagators for B* and 
D* and invariant vertices, and do not rely on an expansion in 1/w,q. In other words, we 
only use the soft-pion Lagrangian to model the pion emission amplitude, assuming that the 
couplings do not vary significantly with energy. In effect, this implies a certain form for the 
contact terms which appear in the chiral Lagrangian at the sub leading 1/ttt.q order. 

In the calculation of the integrated inelastic probability, there is a subtlety that requires 
some care and was discussed already in [10]. Since Mo* > Mr> + m n , the point e = cor- 
responds to \p n \ ~39MeV and the integration extends to small negative e. At e = the 
integral has a singularity related to the Dir decay of the unstable D*, which should be dis- 
tinguished from the actual continuum contribution, and has to be removed. In practice, the 



physical regularization is to introduce the Breit-Wigner factor , replacing 1/e 2 in Eq. (41) 
by l/(e 2 + r 2 /4), where T is the decay width of D*. In actuality T is small compared even 
to the energy release in D* — > Dir. Therefore including the width serves only to regularize 
the integral. Adopting it, integration around e = yields unity, the probability of B^rD* 
we start with, where D* is represented by the Dtt- resonance. The integration over e is then 
carried out with \e\ > e m i n ^> T. The resulting inelastic integral does not depend on the 
choice of e m i n as long as F <C e min <S Mn* — Mn — holds. An accurate analysis shows that 
for all practical purposes the integral can simply be evaluated by setting = Mp* — Mr,. 

It turns out that numerically the most important effect comes from the difference in 
the pion-meson couplings in the charm and beauty sectors, r^l. Various studies suggest 
t i$ 1 [22] • Fig. [7] shows the integrated wP^i as a function of the upper cutoff on the pion 
momentum, p™ ax , for a few values of r. Formally, p™ ax is related to the maximum excitation 
energy e M - 



e M = V(PT X ) 2 + ml + yf (pr x ) 2 + M 2 D - M D * ; 
however, a lower cutoff may effectively be set by the domain of applicability of the soft- 
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Figure 7: Nonresonant Dir contribution including both charged and neutral pion depending on 
the maximal pion momentum, at gn-D-K = 4.9 GeV -1 which corresponds to T D ,+ =96keV. The 
plots show, from bottom to top, r = l, 0.8, 0.6 and 0.4. 
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Figure 8: Nonresonant D*tt contribution depending on the maximal pion momentum, at 
9d*d*it = 4.9 GeV -1 . The solid lines correspond, from bottom to top, to f = l, 0.8 and 0.6, while 
the dashed line refers to f = 1.3. 

pion approximation. It is reasonable to stop at least somewhat below the expected 'radial' 
resonance domain, around p™ ax w 0.6 GeV: even if the amplitude grows with p n , at some 
point this contribution starts to belong to the resonant radial excitations and should be 
excluded to avoid double counting. 

The D*ir channel has been considered only in Ref. [I], even though it is required for 
consistency with the spin-symmetry structure of the corrections already at I/itLq. The 
pattern of the 1/m.Q mass corrections in the intermediate poles is now reversed, and the 
contribution at low pion momentum is suppressed. Altogether it turns out very small unless 
the difference between B*Bir and D*Dn vertices is appreciable. The numerical results are 
shown in Fig. Mrj 

The calculation for D*ti proceeds in a way similar to the Dir case. The 'bare' weak 
5 Due to a typo in the plotting notebook the plot in Ref. [1] showed values larger by a factor of ir. 
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vertices for the mesons are now 



{D\e)\c 1 ^b\B)=2^M B M D , , 
(D*(e)\cw 5 b\B*(e'))=WM B *M D * e„ vpX eV 



Pa 



(42) 



M D * M B * 

where P and P' are the four-momenta of D* and B* mesons, respectively. The pion 



Lagrangian paralleling Eq. ( 39 ) is 



£d*d*tt — 2gD*D*n tfivpx D d u D* d\7c, 9d*d*-k — f 9d*d-k, 



(43) 



from which the D*^-D*tt vertex is derived. The required diagrams mirror Figs. Q, with 
D replaced by D*. The point e = is now below the threshold and no subtleties in the 
integration occurs. Again, the neutral pion contribution to uwi is half of that for the 
charged pion. 

The pion couplings gB*Bir and gD*D*n have not been measured experimentally. There is 
theoretical evidence that l/m c corrections in the coupling should indeed be significant [22J; 
similar conclusions have been reported from the lattice studies [23]. The spin symmetry 
violating difference between Qd*d*-k and go*D-K has not been addressed yet. On general 
grounds it may likewise be thought to be significant; we expect to have its estimate from 
the QCD sum rules. 

We finally arrive at a few percent continuum contribution to 7i ne i: 



'inel 



3 to 5% . 



This is about a fourth of the resonance estimate Eq. (38) and therefore conforms to the 



general expectation. A clarifying comment is in order in this respect. 

The above estimate of the D^n effects should not be added to the principal numeric 



estimate of Eq. (38). The corresponding hadronic states contribute to the nonlocal correla- 



tors along with the resonant states and therefore have been accounted for in the observed 
meson masses used there as inputs. We will come back to this subject in Section 6.3. The 
D^tt states effectively lower the average excitation energy e compared to e ra ,d] however, 
the continuum contribution is relatively suppressed and this is not a prominent effect. On 



the other hand, should one discard the estimate of the inelastic contribution of Eq. (38), it 



is reasonable to include at least in the unitarity upper bound for J-(l); this lowers it 

by about 2% down to approximately J-'(i) < 0.90. 

To conclude the first part of the paper, we refer the reader to Sect. [81 where our numeric 



conclusions for J-"(l) can be found, see Eqs. (131), (132) and (133). 



6 A closer look at the heavy quark excited states 

The zero-recoil transition amplitudes between the B meson and the excited charm states 
appeared in Sec. [5j where they gave the part of the power-suppressed correction to ^(1) as- 
sociated with the overlap of nonrelativistic wavef unctions. They are encountered in various 
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applications of the heavy quark expansion, together with the spectrum of the corresponding 
states, and deserve a dedicated analysis. 

So far mostly the P-wave excitations of the ground-state mesons have been considered 
in the literature. They play a primary role in the small velocity regime and enter the sum 
rules which strongly constrain possible values of the main nonperturbative heavy quark 
parameters [TTI [21]. Here we extend the analysis to higher states, in particular to those 
which are commonly referred to as radial or P-wave excitations, assuming the heavy quark 
(static) limit. In the actual heavy mesons we may have additional excited states associated 
with the interacting spin and gauge field degrees of freedom. 

In our analysis we will use heavy quarks deprived of their spin [TT] [21] ; this is physically 
motivated since in the static limit heavy quark spin degrees of freedom decouple, and this 
greatly simplifies all the expressions, as long as no velocity change is considered. This 
formalism leads in a simple and transparent manner to the trace formalism of the rest 
frame; a discussion of the connection to the conventional formalism can be found in Re. [3]. 
We remind that with spinless heavy quarks the ground-state heavy mesons are spin-| 
particles; the corresponding states are denoted by |f2 ) an d their two-component spinor 
wavefunctions, when needed explicitly, are denoted by \l/o- The Pauli matrices a act on 
the spinor indices of the heavy hadron wavefunctions. To specify our convention we set the 
parity of VLq positive. 

Translating the relations in this formalism to the case of expectation values in actual B 
mesons is straightforward. The spin-singlet quantities are in a one-to-one correspondence. 
For the spin-triplet ones the fact that gq = —a for the spinless state like B meson (<jq is the 
spin matrix acting on the heavy quark indices) immediately gives the required translation. 

6.1 Model-independent spectral representation 

Heavy quark theory requires the tower of heavy quark transition matrix elements with zero 
spatial momentum transfer, ordered according to the number of covariant derivatives (heavy 
quark momentum). They enter both the expansion in velocity and the non-relativistic 1/rriQ 
expansion. In practice we are interested in the matrix elements of the ground state. 

The unit operator QQ has trivial matrix elements at vanishing velocity change. The 
operator Qitj-Q with a single derivative describes the dipole transitions which connect Qq 
to the | and | P-wave states with negative parity. At the lowest order in the expansion we 
need the general inclusive two-point correlator with two derivatives, or its absorptive part: 

-lmP jk (u) = ±- [ d 3 x fdx e^-^(n \Q^ j Q{x)Q^ k Q{Q)\n ) 
it In J J 2M Q 

= T^\uj)^l(8 jk +a jk )^o + T^~\u)*S>l{25 jk -<Tj k )^ Q ; (44) 
here and in what follows we use the notation 
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The two spectral densities generalize the conventional | and | SV amplitudes to continuum 

states: 



UJ) = u 



d|r 1/2 (o;)p 



r<n. 



w = u; 



d|r 3/2 (w)| 2 _ 



(45) 



dw du 
their relation to the inelastic part of the SV structure functions W±(u) introduced in 
Ref. E5| is 



WJuj) 



2T ( t Xoj)+T^ )i 



or 



To extend the analysis to the relevant radially excited states we consider the spectral 
density of the general two-point correlation function of the products of two spatial covariant 
derivatives acting on the heavy quark: 



Rijkl (w) 



1 

2^ 



d 3 x / dx e 



-ILUXQ 



1 



2Mr 



(QalQ^Qix) Q^ k mQ(0)\n ) . (47) 



It can be decomposed into the invariant structures corresponding to three possible classes 
of the intermediate heavy quark states with j = h, j = § and j = §, respectively: 



Rijkl ( w ) 



1 (1+) 1 (1+) 1 (1+) 

-p p 2 [UJ)SijS k i - — ppl {oS){8ij(T k i+<7ij5 k i) + —p g 2 (u)<Tij<Tkl 

y to oo 



1 (§ + ), 



200' 



1 



Pf M ( fiikSjl + $il$jk - + - ( a ik$jl + Vilfijk + CTjkSil + CTjiSik) 



80 Hf 

1 Jl + ), 



16 



(4* 



The meaning of the three invariant functions for j = \ and j ' = | will become clear shortly 
The nonlocal correlators which are relevant to our analysis of Ji ne i have been introduced 
in Eqs. (32). In addition, it is also useful to introduce 

1 



PttG 



d 4 x i\x \ 
d xi\xo\ 



4M B 
1 



2M 



B 



(B\iT{bn 2 b(x), bn 2 b(0)}\B}', 
(B\iT{bn 2 b{x), baBb{0)}\B)', 



\^>s$ijhi + ^P~A(Sik8ji - SuSjk) = j d 4 xi\x \ 



1 



4M B 



(B\iT{baiB k b(x), 6o- i S,6(0)}|S) / . (49) 
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The extra factor i\xq\ compared to Eqs. (32) is simply an extra power of excitation energy 
e in the denominator in the spectral representation, see Eqs. (51) below. Using the p 
correlators one can directly write the inelastic contribution to the sum rule @, cf. Eq. 

j _-(PkG + Pa) ^P-ktv + P-kG , pTTn+PnG + PS + pA { 1 , 2 1 "\ ^ ( 1 

-* inel 



^ )+0 l-T .(50) 



3m 2 2>m c mb 4 \m 2 3m c mb m 2 J lm 3 

All these nonlocal correlators are readily expressed in terms of the above spectral den- 
sities: 

3 f p p 2 \u) _ f pp 2 \u) 

p 3 nn = / du ^ ^ p^ = / dec? 



P^G = / dw ' 7 -P*G = I du 



2pi| +) (^) _ f 2pi| +) (u;) 



3 I a, A Pa' i u ) , Pi' \ - / i, , | Ps 2 M . p/ (wj 



Pi = / do; I =£2L^ _ , p A = / do; 




CJ 2 



,<3 + >,.A ,(f + ), 



2pi 2 V) pi 2 V) 



3tu 2 2u 2 



(51) 



The integration runs over positive u; the point u = is excluded. The integration is also 
cut at large oj for u>> p according to the normalization prescription in the kinetic scheme. 

( - + ) 5 + 

Note that neither Pf 2 j g nor p^ ) can contribute above. Using these relations one can 
express the I/tuq inelastic transition probabilities for actual B mesons. The following form 
appears particularly convenient for the analysis: 



2m c 2m6/ w 2 \2 m c 2mfe/ \3m c to&/ u; 
4 \3m c mJ om~ w z 



The analogous representation for the vector current-induced transitions is given in Eq. ( A. 28 ). 



6.1.1 Intermediate state contributions 

Here we consider the contribution of an individual positive-parity state to the above corre- 
lators. Again, we first briefly remind what happens for the P-wave states. 

Following Ref. [TT] we denote the | and I P-wave states by (p and \ an d describe them 
with the two-dimensional spinor <fi and the non-relativistic Rarita-Schwinger spinors \ji 
respectively. There are successive families of these states which differ by their excitation 
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energy e n = M u — Mq\ we typically omit the index marking the excitation number. The 
|-spinors Xj satisfy o~jXj — and the normalization fixes the sum over their polarizations: 

XiWx]W = Sij - • (53) 

a 6 

The dipole amplitudes are related to the conventional small- velocity amplitudes r by 

(^"Vii^o) = e^f^a^o , (x (m V^o> = v^r^xf (54) 

The two-derivative heavy quark operators at vanishing total spatial momentum acting 
on Vl may create | + , | + or | + states. The ground state itself has j F = \ + \ we will always 
assume however that | + refers to an excited multiplet. We describe the | + and | + states 
by a conventional spinor x an d by the Rarita-Schwinger spinors Xj an d Xjh respectively, 
with the following constraints] 7 " 

a jXj = ; Xjk = Xkj, Xjj = 0, (TjXjk = . (55) 
The sum over polarizations A giving the spin part of the propagator equals to 



$>(A) X t(A) = l 

A=l 

4 x 

^XiWx]W = $ij - n a i a 3 



A=l 



3 1 

/]XijWxiiW = Y^^ikSji+S i i5 jk -l5 ij 5 k i)- — ( (T ik^ji+ (T ii^jk + VjkSii+VjAk)- (56) 

A=l 

Following the standard notation for the diagonal matrix element of the ground state, 

2 2 

(n \n k m\n ) = ^ S kl #J* - ^ *Ja fei * , (57) 
6 o 

we parameterize 

P G 

(| + |7r fc 7rz|^o> = T hi X f *o - 7T xW*o ■ (58) 

The transitions amplitude into | + -states have a symmetric and an antisymmetric structure 
parameterized by constants / and g: 



W = 7^(x\vi + x!^fc)*o + \ itkim xl^o , (59) 



6 The P-wave wave-functions will no longer appear in what follows and we use the same notation for the 
|-hadrons of opposite parity. 
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while the | + amplitude depends on a single parameter h: 

(l + \n k 7r l \n } = h X l i y . (60) 

The residues P,G, f,g and h are different for each multiplet of the excited heavy state. 
A particular hadronic state with excitation energy e n is associated with the following 

factorized contributions to the invariant spectral densities p p f P9>9 (c<;), P// 9g ( w ) an d P^ ( w ) 
in Eq. (pi): 



5p { } + \u) = P 2 5(u-e n ), 5p { p(u) = G 2 5(oj-e n ), 5p p f\uj) = PG5(uj-e n ) 
5pf\u) = f 2 5(u-e n ), 8pf\u) = g 2 5(u-e n ), Sp { p(u) = fgS(u-e n ) 
5 p (I + \uj) = h 2 5(Lo-e n ). (61) 

The ground-state contribution to Rijki is located at u = and is excluded from the nonlocal 
correlators p 3 and p of Eqs. (51); it is obtained by taking P( ) = p 2 and G( ) = Pq with 
eo = 0. 

The spin- 1 contribution 5p^ \u) is given by the square of a single residue h and is 
non-negative. The factorization relations lead to the Cauchy inequalities 

(pifV < pf ] pf\ (P/fV ^ pPpP (62) 

which turn into equalities for a particular single state contribution. Furthermore, one has 
the general inequalities 

2\p { P f\u)\ < pp +p f\ 2\pf g + \u)\ < pf ] +pf\ (63) 

the first of which will soon be useful. 

For completeness we mention the constraints imposed by the BPS limit for the ground 
state, which assumes na\Q ) = 0. Similar to the ground-state expectation values, for transi- 
tions into ^ + states the BPS condition implies P — G\ however the relations for the excited 
states are accurate only to first order in the deviation from the BPS limit. No constraint 
emerges on h, while for the transitions into | states BPS implies f = g: 

(l + \n k na\n ) = (lf-lg)x{y - (64) 



For the invariant structures in the spectral density Rijki(to) in Eq. (48) the BPS condition 
leads to the relations 

pf\co) = p p f\u) = pf\u), pf\u) = pf\u) = pf\u). (65) 

Let us note that when the spin-^ degrees of freedom associated with light antiquark 
in the meson decouple, as in the case of purely perturbative calculations, the correlator 
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decomposition can be based on the angular momentum only, L = 0, 1, 2, of the intermediate 
states composed of the QCD degrees of freedom still interacting with heavy quark: 



9 QijOfei H {OilOjk-OikOjl) H j^(OiJfe<5j7 + <5«Ojji.--<Jy<3fciJ 



^o, (66) 



with a trivial spinor structure. This would mean 

Pp \u) = PoM, V) = f pi(u;), pg 5 = px(u), 



pf\u) 



iP2{w), p ( t +) {u) = ^p 2 {UJ), r/(/ 
We will use these relations when address the perturbative p-dependence later on. 



i 



j§ + ), 



Pfs [UJ] 



0. 



(67) 



6.2 Hyperfine splitting and the estimate of the correlators 

The analysis of the spin structure and of the factorization properties does not constrain the 
overall scale of the residues and the significance of the radially excited states, which can 
however be estimated by considering the heavy quark mass dependence of the hyperfine 
splitting AM 2 . This will allow us to fix the magnitude of the nonlocal correlators in 



Eqs. (31), (50). 



The B — B* mass splitting basically fixes the value of p? G . The mass difference between 
the beauty and charm states mainly reflects nit, — m c . A comparison of Mp* — Mo and 
Mb*— Mb hyperfine splittings gives information on higher dimensional correlators, and in 
particular on the values of two D = 3 spin-triplet parameters, p\ s and p^ G + p\: 

2 Pg , 'Pis + pIg + Pa , 1 



AM B = M b -Mb = -^+ - r f ' +Ol—A (68) 

6m b 6ml \ m b J 

and likewise for charm. Following Ref. [6], we explore this relation cast in a somewhat 
different form. 

As dictated by the heavy quark expansion, for sufficiently heavy quarks the difference 
AMq = Mq,—Mq for the vector and the pseudoscalar mesons approaches a constant. Yet 
it has been noticed long ago that such a relation works well even for light quarks: 

M 2 -M 2 ~ M^-Ml ~ M 2 D *-M 2 D ~ Af|. -M 2 B . (69) 

Clearly, the heavy quark expansion cannot explain why such a relation extends down to the 
light quarks. The approximate equality resembles the universality of the slope for Regge 
trajectories and may root in the peculiarities of the strong dynamics. Moreover, actual 
QCD predicts that such a relation must be violated for sufficiently heavy quarks due to the 
perturbative renormalization of the chromomagnetic operator of the heavy quark which has 
a nontrivial anomalous dimension. As was also noted long ago, the observed 12% decrease 
in the mass square difference in the B system compared to D mesons fits reasonably well 
the naive estimates of this perturbative evolution. 
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It is therefore tempting to assume the independence of the mass-square splitting as a 
phenomenological property of soft strong dynamics yet to be understood, and to consider 
its consequences in the context of the heavy quark expansion where we can vary the heavy 
quark mass. Using the mass formulae we have 



AM 2 =AM Q (2M Q+ AM Q ) = \c G {^ m Q )& + \ p1g + P " A Pls + + 0\ '. ). 70 i 

6 6 mq 




where cq is the Wilson coefficient of the chromomagnetic operator in the heavy quark 
Hamiltonian. If the perturbative mg-evolution of cq accounted completely for the observed 
difference between AM^ and AM|>, the m^-independence of AMq to the first nontrivial 
order in l/mg would imply 

- {~pIs + pIg + Pa) - 2A^g • (71) 

This relation will guide our analysis below. 

Numerically using A = 600MeV, /4 = 0.3 GeV 2 and p 3 LS = -0.1 GeV 3 we would get 

- (PIg + Pa) « 0.45 GeV 3 (72) 

which indicates that the nonlocal correlators p 3 are in general sizable. The negative sign 
for pl G +p\ is expected from BPS arguments [6]. 



To proceed more quantitatively we introduce a factor k in Eq. (71), to account for the 
actual mismatch between the observed mass dependence and the dependence stemming 
from cq\ 

~ {-As + PIg + Pa) = 2(l+«) Afi 2 G ■ (73) 

k does not include higher power corrections, which we address later, and can be defined 
through 

~3 



nAp G = lim m 2 



A AM*-c G (m Q )p 2 G 



(74) 



where the logarithmic derivative of cq^g ls related to the exact anomalous dimension of 
the chromomagnetic operator, times AMq. 

In simple words, the soft part of AMq, identified by subtracting the (Wilsonian) short- 
distance renormalization factor, approaches a finite limit as mq — > oo. The assumption 
underlying the approximation of small \k\ is that such a soft part is nearly a constant in a 



wide range of heavy quark masses, as Eq. (69) would naively suggest. 



Let us now look at the perturbative renormalization cq- The one- and two-loop [26J 
contributions are known. At first glance in the evolution from beauty to charm, namely 
in the ratio CG(rn c )/cG(rn b ), the two-loop contribution is quite sizable. However, the bulk 
of the higher-loop perturbative enhancement comes from growing strong coupling at small 
momenta of virtual gluons. This large-coupling domain must be removed from the pertur- 
bative corrections to avoid double counting. The subtraction piece is power-suppressed yet 
important for charm, especially in the effect of running of the strong coupling. 
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CG(m Q )/c G (m b ) 




1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 



Figure 9: The effect of removing the infrared piece belonging to the nonperturbative expectation 
values, from the one-loop perturbative evolution from beauty to charm CG(mQ)/cG(rrib), as a 
functon of tuq. The solid curves correspond to \i = 0.8 GeV (blue), 0.7 GeV (red) and 0.6 GeV 
(green), respectively; the dashed curve shows the naive perturbative result with no cutoff, all at 
a s = 0.3. For comparison the dashed-dotted line shows the no-cutoff result at a s = 0.22. 



The numeric aspects are illustrated in Fig. [9] in the case of one-loop corrections, using 
a s = 0.3. The cutoff effects are more important than higher-loop corrections. The two-loop 
calculation of Ref. [26], which strongly enhanced the first-order renormalization, should 
not be used in the present context To get an accurate estimate in the following we 
will employ the one-loop calculation of cg{tti c ) / CG{jn b ) with Wilsonian cutoff evaluated at 
a s = 0.3. 

A precise determination of k is not easy since only two data points exist on the AM 2 
curve in the heavy quark regime. Moreover, the perturbative treatment of charm may have 
insufficient precision and higher power corrections for may be significant. The situation 



is illustrated in Fig. 10 where only the change in the perturbative coefficient cq is consid- 



ered. In other words, the continuous lines show the quantity (Mg* — M|) i without 

power corrections associated with using the meson masses Q Clearly, this approximation is 
valid for sufficiently large masses. 



Let us now discuss the implications for the nonlocal correlators that follow from Fig. 10 



The solid curves correspond to k — 0, and do not include any power correction. For K^Oa 



power correction appears according to Eqs. ( 70 ) on top of the perturbative renormalization. 



Therefore, if the actual AM 2 below beauty goes lower than the perturbatively continued 
dependence, k is positive; conversely, a negative « corresponds to the case where the actual 
AM 2 increases steeper than the one computed perturbatively from the beauty point. We 
also note that what matters here is actually the relative position of the curves in the vicinity 
of the beauty point, or, more generally the difference in the slope of the two curves at any 



7 The physical masses depend on the pseudoscalar meson mass Mp, viz. Ms, Mp>, ■■■ whereas the per- 
turbative renormalization is expressed through the quark masses. To put them on the same plot and to 
compare we evaluate the perturbative calculations at tuq = Mp Q — A. The dependence on the precise value 
of tuq is minor. 
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Figure 10: The difference AM 2 = My — Mp plotted as a function of the inverse meson 
mass 1/Mp. The three data points stand for beauty, charm and strangeness, the horizontal 
line shows the asymptotic value given by the p-ir splitting. (The K point has been moved left- 
ward.) Solid curves represent the expected perturbative evolution in the case all power correc- 
tions are neglected, using the Wilsonian perturbative renormalization factor; they correspond to 
\x = 0.8 GeV, 0.7 GeV, 0.6 GeV, respectively. The dashed curves shows the naive perturbative 
result without cutoff. 



sufficiently large mass where the 1/m.g and higher terms in the masses can be neglected. 

It looks improbable that the actual AM 2 dependence on the inverse meson mass may 
be significantly steeper around the B meson than the one-loop a s = 0.22 perturbative curve, 
because that would require an unnatural shape, perhaps with a maximum higher than the 
charm point. Moreover, in that case sizable power corrections would be necessary to hit 
the AM|, point, and there would be no reason to expect a small value of the nonlocal 



correlators. Since the AM|, point in Fig. 10 is above the Wilsonian perturbative curves, a 



somewhat steeper dependence on 1/Mp is observed. Therefore, n must be relatively small 
and negative in the Wilsonian approach. 

In order to obtain a numerical prediction for AM 2 and refine Eq. (72) we use the ex- 



pansion for AM, Eq. (68), and rewrite it, neglecting, as usual, the perturbative corrections 
to power corrections, as 



m b c G (m Q ;fi) 



5_ 

a !»■/ x ni b c G yuiQ]H) - mQ 

AMimQ ^ mQ ca^)i-^ AM{mb n 



O 



i 

mi 



(75) 



with the shortcut 



8 



-(pIg + Pa-Pls) 
2/4 



or 1 + k 



5_ 
1 



Higher order power terms modify the mg-dependence of AM 2 . Therefore, in order 
to obtain definite numerical values for the nonlocal correlators from AM 2 in beauty and 
charm we need to make assumptions on the higher-order terms Oil/rnq) in Eq. (|75j). The 
simplest option is to discard 1/mg and higher terms in the mass difference. An alternative 
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Figure 11: AM 2 for different heavy quark masses computed under different assumptions for 



higher-order power terms: no higher power term (left) and the ansatz in Eq. (76) (right). The 
two families of curves assume different values for 1 + k from top to bottom: 0.5, 0.55, 0.6 (left) 
and 0.7, 0.8, 0.9 (right), respectively. A = 0.65GeV, /x = 0.8GeV and a s = 0.3. 



way is to write the power correction factor as 



AM(m Q ) = 2 -c G (m Q )^.-^ r (76) 

O Tfln 1 H 

as an ansatz for higher-order terms; it expresses them in terms of /j/q and of the powers 
of 5. The two forms have identical /^h a dr/ m Q corrections but differ at higher orders. The 
latter ansatz has the advantage of yielding a reasonable finite value even at small TOq. 

To depict results graphically as a function of Mp, for the pseudoscalar mass Mp^tjiq) 
we use its conventional heavy quark expansion well described by just rriQ + A due to the 
proximity to the BPS regime (for even lighter quark, to continue the plots below charm we 
use an ad-hoc extrapolation giving a reasonable behavior at small mQ). 



Fig. |11| shows the expected AM 2 for three values of « under the two assumptions 
about higher-order terms. The difference illustrates the sensitivity to the unknown power 
corrections in the meson masses. It is modest around the B mass scale, yet becomes 
significant for charm. We have also tried alternative ansatze for power corrections and found 
that k always varies between —0.5 and 0. The extrapolation of the curves below charm 



seems to suggest a preference for the smaller \k\, as in the right of Figs. 11 Moreover, higher 



order perturbative effects are qualitatively expected to increase 1 + k. We will therefore 
use k = —0.2 ± 0.2, arriving at 

2(1 + k)Ah% ps 2(0.8 ± 0.2) 0.65 • 0.35 GeV 3 ~ (0.35 ± 0.1) GeV 3 . (77) 

Combining this estimate with the expected value p\ s ~ —0.1 GeV 3 we end up with 

-(p'c + Pi)^0.45GeV 3 , (78) 



as in the original estimate, Eq. (72). One of the reasons is that a larger value for A must 



be used in the Wilsonian approach, which offsets a negative k. 
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The main assumption employed in the above estimate is the possibility to quantitatively 
use the mass expansion for charm particles, assuming a reasonable magnitude of the higher- 
order power corrections in a physical scheme. This assumption could be avoided in the case 
an additional input from lattice calculations were available, with the squared mass splitting 
reliably computed at one or more intermediate points for quarks heavier than charm. The 
short-distance expansion there is in a better shape, perturbative corrections under better 
control, and higher-order power corrections are less important. This would allow to use 
the differential version of the constraint on the hyperfine splitting in the large-m^ limit, 



namely Eq.(74), fixing directly the slope. Even if a first-principle direct measurement of 
AM 2 with sufficient precision turns out difficult for large mass, any intermediate point - 
and even a point somewhat below charm - would constrain the shape of the curve thereby 
narrowing down the interval of possible values for the derivative. 

Eq. ( J78| is a strong constraint, especially useful in the analysis of the inelastic contri- 
bution to ^"(l). This will be discussed in the next subsection. Here we only note that a 
determination of — (p^-g + Pa) gi ves a model-independent constraint on the spectral densi- 

(- + ) 

ties pg , ppl and p g 2 , which are informative due to the factorization property and to 



the relation to the spectrum of physical states. Since the value in Eq. (78) is quite large, 

p P g (oo) together with — p 3 G is most probably positive, unless the correlators are strongly 
dominated by the | + states. The predictive power sharpens significantly if additional dy- 
namic information is used, for instance the proximity to the BPS limit or the saturation by 
a limited number of the heavy quark multiplets. 

We conclude this section noting that another interesting constraint on the nonlocal 
correlators can be obtained from the difference of the spin-averaged B and D meson masses, 
namely from 

m b -m c = M B -M D + & (—-—) + ( ~ 2 ~ ~ 2 ) + O [ — ) , (79) 

where p 3 = p^. 7r +p|. The use of this relation to determine p 3 requires a good control of the 
heavy quark mass difference and a precise knowledge of p 2 , and p|,, which can in principle be 
provided by the global fits to semileptonic moments. We have considered the fits employed 
in Sec. El they typically give p 3 = (0.33 ± 0.17) GeV 3 . Notice that the semileptonic fits 




determine p 2 in the actual B meson, while its static approximation appears in Eq. (79). As 
the difference is BPS and 1/rrib suppressed, it is certainly smaller than the fit uncertainty 
on pi, and we have neglected it. The value of p 3 can be linked to (p 3 G + p A ) since the 
sum p 3 ^ + p| + p 3 G + p\ is positive definite, vanishes in the BPS limit, and receives only 
second order corrections to the limit. Therefore, — (pI-c+Pa) ~ (0.33±0. 17) GeV 3 , which is 



compatible with our primary estimate, Eq. (78). 
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6.3 The hyperfine constraint and the excited states 



The hyperfine splitting constraint Eq. (78) fixes the overall scale of the correlators which 
are the focus of our study: 



do; 



0.45 GeV 3 + (k+0.2) ■ 0.36 GeV 3 . (80) 



The normalization point enters here as the upper cutoff in the integral over u; we assume 
it is around 0.8 GeV. At first glance we have a single relation among four spectral func- 
tions; nevertheless this relation turns out to be very useful thanks to the positivity and 
factorization properties discussed in Sec. 6.1.1. 



The principal observation is that according to the estimate Eq. (78) the nonlocal corre- 



lators are numerically large. The particular combination constrained significantly exceeds 
the estimated size of the expectation values of local heavy quark operators of D = 6, most 



notably p D which by itself is large. Alternatively, the value in Eq. (78) can be compared 
to A : it exceeds it, even though A is already numerically large. Likewise, the square of 
the kinetic expectation value (p^) 2 can be used to gauge the scale of hadronic parameters 
of mass dimension 3 if it is divided by a typical energy of the radial excitations, yielding 
again a value close to 0.35 GeV 3 . 

For more quantitative statements it is advantageous to separate the resonant and the 
soft continuum contributions to the spectral densities: the inclusion of the continuum makes 
the effective value of u identified with e ra d more uncertain numerically. The separation can 
be done since the Qq+tt states are readily analyzed in the soft-pion technique; for instance, 



Sect. 5.2 estimated the continuum part of 7 ine i at actual values of m c and m&. This analysis 



is performed in Sect. 6.4 and suggests that a relatively small fraction should be subtracted 



from the value in Eq. (80) if we want to keep only the resonant contribution. 

Once only the resonant contribution are retained, a meaningful approximation well 
illustrating the physics of the constraints is a model where a single multiplet in each channel 
involved is considered. This model has a very few parameters essentially reduced to the 
residues P, G and g (/ and h do not enter), and all the constraints can be easily analyzed. 

To gauge the size of the related hadronic parameters we take the sum of the two spin- 
singlet nonlocal correlators p^ n and p| as a measure of their overall significance; each of 
these correlators is positive. On the other hand, we have a relation 



p 3 



+ Ps 



-(pIg + Pa) 



doo 

CO 



Pp 



2p 



P9 



Pi 



(81) 



The last bracket on the right is positive; it can be represented as a correlator of a certain 



operator. This shows that the l.h.s. of (81) always exceeds the combination —p 3 T c~PA fixed 



by the 'hyperfine' constraint; the minimum is attained at the BPS point where the two 
combinations coincide. The inequality actually holds at arbitrary spectral parameter u; it 
is the numeric hyperfine constraint that applies only upon integration over the energy. 
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Eq. (81) shows that the 1/mg spin-averaged meson mass shift must always be larger 
than 3/4 times the I/tuq correction to the hyperfine splitting. This is not difficult to trace 
directly: 

"«l+tHB = * {iM ._ iMB)+iMB . (82) 

5i/ m ^MB receives a positive contribution (Pd+Pls)/^ 171 ! fr° m the local piece and a negative 
contribution from the T-product of two (<jtc) 2 . Both vanish at the BPS point. 



The last positive term in Eq. (81) is of the second order in the deviation from BPS, 
therefore it may be a good approximation to neglect it in numerical estimates unless the 
BPS is strongly violated in the actual B mesons. For what follows we do not need this 
additional assumption, the positivity is sufficient. 

Now we can see that the sum of the two correlators p\ % and p\ is also quite large 
numerically; this implies large 1 / uiq corrections to the meson states, in particular averaged 
over the spin multiplet. However, experiment tells us that the heavy quark symmetry works 
reasonably well even in charm. This means the correlators driving power corrections should 
not be excessively large. In order to have both p^ n and p| as small as possible, besides the 

1+ 3 + 

BPS condition, we can assume that pl n and p| are equal, in which case p| = \p} must 
hold upon corresponding integration over energy. 



On the other hand, since the r.h.s. of Eq. (80) is large, it is reasonable to expect that 



p p 2 g is positive, and that p p 2 is not much smaller than p g 2 . In this case the left-hand 
side is a sum of two positive contributions from ^ + and | + channels. 

It makes sense to parameterize the relative contribution of the ^ and | states, intro- 
ducing a parameter v. 

du pf\u) _ fdu ( pf\u) (|+) 2pf\u) 

— g — / ' — ^ — + p p 9 ( ' * — I ' ^ ' 



The large numeric value of the the r.h.s. of Eq. (80) points at v between and 1 as the 



most natural solution; v exceeding unity is highly improbable. 
6.3.1 Estimate of I x 

The spectral representation together with the hyperfine constraint allows one to analyze 
the possible values of L me \. Ignoring a possible spread in the values of average excitation 



energy we equate / ine i with Ii/e ia a or, equivalently, express the p correlators in Eq. (51) 
through the conventional p 3 . The evaluation of 7i ne i then depends directly on the size of 



the corresponding nonlocal correlators. For instance, the last positive term in Eq. (31) is 



immediately recognized as the square bracket in Eq. (81); it is given solely by the BPS- 
violating transitions into the | multiplet. 

All the terms in l\ but the leading BPS piece fixed b y th e hyperfine constraint are 

o _|_ 

independent of the § contributions, and the l.h.s. of Eq. (80) gets a positive piece from 
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them. Therefore the minimum of I\ is attained at vanishing p g 2 which also corresponds 
to p\ = 2pf . 



We now refer to Eq. (52 ) for the structure of Ji ne i and express it in terms of the transition 
amplitudes P, G, g. At fixed G and g, Ii ne \ is a quadratic trinomial in the ratio Pj G, with 
the positive leading power coefficient proportional to (l/m c — 1/mf,) 2 . We recall that in the 
BPS limit P/G = 1. The hyperfine splitting implies an additional constraint on P, G, g, 
which can be taken into account. The result is that 

:i '"' l> > 1 - (1 - z/)^f ~ 1 - 0.07(1 - !/), (84) 



-(pIg+Pa) 
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and I\ has a minimum at a value of P/G > 1, whose exact position depends only on 
the heavy quark mass ratio but not on g. The minimum is attained, with the inequality 
saturated, at 

Goo m^—m c 

wher e th e numeric values correspond to m c /mf, = 1.2/4.6. This result was referred to in 
Sect. O For v = it implies G « 0.37 GeV 2 . 



The minimum, however, is rather shallow, see Fig. |5j because the coefficient of the 
term quadratic in P/G is large while the linear term is small. While I\ may, in principle, 



significantly exceed its BPS value If PS , Eq. (36 ), for strong BPS violation if P/G is negative, 
at more natural positive Pj G (a regime where both p\^ and p\ are not too large) the value 
of Iijlf PS is typically slightly smaller than unity, as illustrated in Fig. [5] 

1 + 3 + 
6.3.2 The inclusive ■= and 5 yield 

The hyperfine constraint implies relatively large values of the transition matrix elements to 
| + and/or | + states and an enhancement of the overall yield of the corresponding excited 
charm states in the semileptonic decays of actual B mesons. This is quantified by the 
estimate of / ine i constituting, roughly speaking, 15% of the D* probability, in the zero- 
recoil kinematics. This comparison does not include the phase-space suppression of excited 
states resonances which becomes significant when one integrates over all available phase 
space. Other kinematic effects may work in the opposite direction, however, and a more 
substantiated estimate is desirable. 

Our preceding analysis constrains directly the 1 / rriQ squared transition amplitudes into 
the corresponding hadronic states at zero velocity transfer. At first glance, this seems too 
crude an estimate at non-zero recoil, where contributions not suppressed by 1/m.Q appear. 
However, since we consider states with the quantum numbers of radials or of -D-waves, the 
leading-order amplitude is proportional to the second power of the velocity of the charmed 
meson — roughly speaking, it is generated by replacing the momentum operators iTk with 
rriQVk- Since for excited mesons Vk is always relatively small, the zero-recoil amplitude 
receives only a minor correction when integrated over the whole phase space. Therefore, 
we can retain only the 1/mg part of the transition amplitude. For comparison, the leading 
heavy-quark transition amplitude into the P-waves provides only a suppressed correction 
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to the 1/ttlq term [27J, even though the leading amplitude in that case is of the first order 
in the velocity, and the P-wave states are lighter. 

We will therefore relate the transition probabilities into radially excited states directly 
to the zero-recoil observables we have studied in the previous sections. The total yield is 
also fed by vector-current transitions, even though we do not expect them to provide a 
large contribution, as their amplitudes vanish in the BPS limit. These transitions can be 
considered in the same way as the axial timelike component, see Appendix |D} 

The total yields are roughly proportional to (Mb — M) 5 = A 5 where M is the mass 
of the corresponding charm multiplet. Numerically we equate A with M B — M D — e ra( j, 
leading to A ~ 2.6 GeV. However, for A/Mb ^0.5 relativistic effects modify significantly the 
nonrelativistic A 5 dependence. A more accurate approach is the following (see Appendix |d| 
for details). Instead of considering the transition amplitudes between the B meson and 
all the excited mesons belonging to a multiplet, we evaluate the decay rate of the VLq 
spin-^ heavy state into the corresponding excited half- integer spin multiplets. The weak 
current coupling of these fictitious hadrons is then fixed by the corresponding transition 
probabilities near zero recoil. 

In practice, we parameterize the transition amplitudes to | and | states, for axial and 
vector currents, and compute the decay rate integrating over all the available phase space. 
The overall normalization of the amplitudes is fixed at zero recoil by the corresponding 
nonlocal correlator. Let us illustrate this in the case of the transitions into ^ + states. 

The most general vector or axial vertices in terms of full Lorentz spinor wavefunctions 
have the form 



M B '° "M B 



o- 



J V , = 9vX^% + bvX^^o + c V j^X%- (86) 

Mb m b 

At zero recoil the effect of h A reduces to a change in qa of (Mb—M)/Mb &a, while the term 
proportional to ca vanishes. However, for generic recoil these are all independent structures. 
Similar considerations apply to the vector current. We neglect the additional contributions 
at non-zero recoil and effectively set the formfactors bA,v, ca,v to zero. Likewise we neglect 
the velocity-dependence of the formfactors. The corresponding decay rates are then given 
by 

1 v,A - \9v,a\ ig2T^ z A,v{r), (87 J 



where 



z A ,v(r) = z ^^~ z ^ r \ 2 (r) = l-8r + 8r 3 -r 4 -12r 2 lnr, (88) 



2 

M 2 



M 2 ' 



zq{t) = 2y/r [l + 9r — 9r 2 — r 3 + 6r(l+r) lnr] 



are the weighted phase space factors for the axial and vector transitions, respectively; z (r) 
is the standard kinematic factor for V — A decays. 
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The value of \qa\ 2 is directly given by the | contribution in w- me \{u) in Eq. (52) for 



one multiplet, and the transition amplitudes G and P are constrained by our analysis of 
the hyperfine splitting. \gv\ 2 follows from the analogous relation for vector- induced prob- 
abilities, Eq. (A. 28). The vector-induced probability is significantly suppressed compared 
to the axial one. 

As already mentioned, we neglect additional recoil corrections. A justification for this 
is the 'extended' SV regime relevant in the context of a large-5 expansion jl] where the 
enhanced corrections reside in the lepton phase space. The deviations from the small-recoil 
kinematics in the amplitude are not enhanced, but rather suppressed by the large power of 
energy release decreasing the average recoil. 

The case of the decays into the | + is treated similarly — the details are given in Ap- 
pendix |Dj We describe these states by complete relativistic Rarita-Schwinger wavefunctions 
at arbitrary velocity, and calculate the corresponding contribution to the (unpolarized) 
structure functions of Qq. Their integration yields the total decay rate, and we fix the 
normalization of the formfactors at zero recoil. 

In this way the decay rates are calculate separately for the | + and | + states and sepa- 
rately for the axial and the vector transitions. Numerically we get for the combined yield 

r 

0.07, (89) 



r s i 

with the axial part strongly dominating (we have included perturbative corrections in the 
denominator). While the relative weight of | and | states varies depending on their 



couplings, their sum is approximately fixed by the hyperfine condition (80). 

In our approach the total decay rate is computed in the heavy quark limit. We in- 
cluded the 1/rriQ symmetry-breaking terms which mediate the transitions in question, and 
exploited the fact that the total decay rate does not depend on the exact mixing of the 
final heavy-quark eigenstates in the actual QCD hadrons. 

To apply the estimates to the yield of actual individual charmed mesons in QCD one 
would need to properly construct the corresponding states in terms of | and | states. 
While straightforward in practice, this may yield unreliable predictions in practice since 
the heavy quark symmetry is probably quite strongly violated for the excited charm states. 
We expect much more robust predictions for the total yield summing up all the associated 
channels. Such inclusive probabilities are insensitive to the details of the strong Hamiltonian 
in the final states, are not affected by possible degeneracies and altogether enjoy smaller 



preasymptotic corrections. Therefore, we view Eq. (89) as a good starting estimate of the 
overall yield of the descendants of the 'radial' states. It refines the earlier estimate given in 
Ref. P], and probably represents a natural lower limit. A compatible number has recently 
been suggested in Ref. [28J. 

It is worth noting that the semileptonic phase space factor strongly suppresses the yield 
of the states with higher mass; in particular, for a wide resonance with mass M = M+5m 
one has 

5 m 



(M B -M) 5 = (A-Smf ^ A 5 1 + ... (90) 

1 600 MeV 1 



37 



This shows that for broad resonances the phase space factor averaged over the whole de- 
cay kinematics significantly distorts the Breit-Wigner shape and shifts the apparent peak 
towards lower mass. The factor in Eq. (90) applied to a wide resonance may easily mimic 
a typical non-resonant continuum yield with a threshold suppression. 

The above estimate suggests that the total yield of the discussed 'radial' excitations 
with mass below M D + lGeV is expected to be, in terms of T sl , at the 7% level. This 
is close to what is observed in experiment, yet traditionally is attributed to the 'wide' | 
P-wave states. The conventional allocation creates a problem: theory predicts that the 
|-states must strongly dominate among the P-waves; they have been measured at the right 
rate of about 10% of T sl , see e.g. |29j for a discussion. We are therefore led to argue that 
the decays into the | states are indeed suppressed, while the bulk of the experimentally 
observed "wide" structure is actually the result of the significant fraction of the 'radial' 
states. 



6.3.3 The radial and D-wave excited states 

The phenomenological analysis earlier in this section suggests, as the most natural solution 

1 ~\~ "f" 

that both the | true radial excitations and the | P-wave states contribute significantly 
the spectral functions involved and are produced at appreciable rate in the semileptonic B 
decays. At present we do not have accurate enough data to state which of the two channels 
dominate; we have only observed that the solutions where one of them, say, the P-wave, is 
small is disfavored. 

In fact, the identification of the two families of states with the S- or P-waves is not 
strict in the actual mesons in QCD even in the limit of a large mass tuq. For instance, the 
| + states can instead be a result of excitation of gauge degrees of freedom. Whether this 
is so or not is an open question, and the fact that we expect a numerically large transition 
matrix element of the chromomagnetic field, 

(l + \QB k Q\n ) (f= o = ^xl^o, (91) 

is not directly related to it. It appears that so far the radial states attracted more attention 
in the literature, while the P-waves were marginally considered [30J . At the same time we 
expect the dominant mechanism for the production of the corresponding charmed states to 
be the l/m c -component of the amplitude, and it is qualitatively similar for both of them. 
Their masses are also expected to be in the same range about 700 MeV above the ground 
state. 

A challenging question is how one can disentangle these states in experiment. In the 
simplest constituent quark model one expects the hyperfine splitting inside the P-wave 
multiplet to be particularly suppressed. However, it is not evident to which extent this 
would be the property of the actual QCD states. Moreover, the hyperfine splitting may 
well happen to be suppressed within the radially excited multiplet as well. 

Some differences are expected in the decay pattern. (We reason in terms of the asymp- 
totic states deprived of the heavy quark spin; the translation into the actual mesons is 
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standard.) The radially excited states can decay into the ground state Qq and a single pion 
in the P-wave, or into Qq plus two pions in S-wave. We expect the dominant channel to be 
the latter where two pions have a cx-meson enhancement; that is, they must predominantly 
be in the isospin-singlet S'-wave state. 

This particular two-pion channel is not allowed for the | + state which should then 
decay mostly into Qq and a pion. A weaker two-pion channel cannot show the resonance 
enhancement associated with a-meson. We naturally expect the | states to have a smaller 
width. It is possible that the excited mesons recently reported by BaBar [31 J with mass 
around 2.75 GeV are related to these states. The states with the lower mass around 2.6 GeV 
may be the radial states. 

Another decay chain where the first decay proceeds into a single pion in the S'-wave and 
a P-wave charm state, either | or | , can be competitive and may provide an additional 
handle through the identification of the P-wave state via its subsequent decay. These 
questions deserve further consideration. Since the two multiplets are expected to have close 
masses, the actual excited vector D* states may show significant mixture, which has to be 
considered. 

6.4 Nonresonant D^iz in the spectral representation 

The special role of the nonresonant D^tt states manifested itself already in the analysis 
of J-"(l). In the heavy quark limit we consider their counterpart, the non-resonant states 
flo + 7c (kaon or r] may also be included) . They are of special interest for a relatively soft 
pion where its energy is essentially below the resonance excitation gap. This gap depends 
on the orbital momentum of the pion: for the P-wave states it is about 400 MeV. Our focus 
is on the radial (or P-wave) excitations where it is about 700 MeV. 

The traditional classification over spin-parity of the light degrees of freedom is equally 
applicable to multi-particle states, including the VIqtt continuum. They can be classified 
in a way similar to the ground-state excitations, the analogies of P- waves etc., and only 
have to be additionally labeled by the continuum excitation energy. In the static limit the 
latter is equal to the pion energy. The quantum numbers of QqTT are not fixed a priori: a 
continuum state is generally a mixture with different quantum numbers depending on the 
production amplitude. For instance, this would apply to the relative weight of the | and 
| states in P-wave. 

The advantage of the expansion we employ is that only the heavy quark states with 
vanishing total spatial momentum are involved, and they are considered in the static limit. 
This fixes the structure of the transition amplitudes appearing to a particular order in the 
1/rriQ expansion; the relative weight of different spin-orbit multiplets is then determined as 
well. In this section we obtain the decomposition into |, | etc. states for Qqtt and calculate 
the corresponding spectral densities as a function of the pion energy. They then determine 
the contributions to | t~i/2 1 2 , l r 3/2| 2 , pL- e ^ c - F° r instance, we will see that the pion loop has 
a T "i/2 = r 3/2 property [32] ■ The associated continuum states do not contribute to fi G or p 3 LS , 
yet they change the IW slope or A in a predictable way, and mediate a positive contribution 
to 
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Figure 12: The pion emission diagrams in the static limit. The solid blocks denote the QiDjQ 
or QiDjiDiQ operators with q = 0. Pion carries spatial momentum k; incoming 0,q have vanishing 
spatial momentum and outgoing —k. 

The Qqtt states at rest are uniquely characterized by the energy and by the pion orbital 
momentum L. Indeed, the total angular momentum j consists of | of f2 and of L of pion: 
L = j ± ~. Its parity relative to parity of is (— 1) L+1 . Therefore, the combination of j 
and parity unambiguously specifies L. For instance, L = are | P-wave states, L = l give 
^ + and | + 'radial' excitations and the § P-wave states require L = 2. 

The combination of the 'P-waves' appearing to the leading order in either l/rng or in 
velocity is mediated explicitly by the operator QiDQ; the relative mixture is determined 
by a concrete form of the amplitude which involves the spin of Q®. The most straightfor- 
ward approach is to consider the Qoir contribution to the zero-recoil correlation function of 



operators QiDjQ and QiDkQ, as we did in Sect. 6.1, Eq. (44). 



The effective low-energy Lagrangian of the 7rf2 f2o-interaction corresponding to Eqs. (39) 



and (43) in relativistic notations is 

£ c hi = -gB*Bir&o7 t i75tt d tl 'K = -2g B * BlT M Q Q, i'y 5 tt n (92) 



(total derivatives are omitted). The diagrams to be calculated are shown in Fig. 12 where 
the heavy hadron lines now all refer to Qq and the solid vertex stands for the operator 
QiDjQ. The vertex is simple: 

' (^o{P2)\Q{iDj ±iD j )Q\n { Pl )) = {p 2 ±p 1 ) j ^ , (93) 



2M Qo 

where we have generally distinguished the left and right derivatives for the case of different 
momenta of Qq- Since in our case the spatial momentum flowing into the vertex vanishes, 
this specification is superfluous. 

The diagram a) vanishes, while b) yields an amplitude with a simple spin structure: 

' ■(n Q *\QiD j Q\n (p)) = 9b* Bw - tfjid^o, (94) 



2M Qo u 

where k is the pion momentum. The resulting correlator Pj^iuj) obtained by squaring the 
amplitude and summing over polarizations of intermediate Qq has only a part symmetric 
in j, k: 

-lmP jh {u) = %^%5 jk ^ Q . (95) 

7T lZ7T z UJ Z 



40 



From this we read off 



T^\u) = T^\u) = %*f6{uj-m^-mlfl* (l-^) , (96) 
where the overall factor refers to a single charged pion loop contribution. The neutral pion 



additionally contributes a half of that. Relations (45) give the values of the corresponding 
\ti/2\ and |r 3 / 2 |, which are equal. 

It is worth noting that the equality ti/ 2 = t 3 / 2 is not an automatic property of quantum 
numbers in the Dtt system. It rather follows from the form of the soft-pion amplitude for 
heavy mesons. For instance, a structure ^oCj^o would have produced only | but not | . 
The relation generally changes already due to the final state interaction (FSI) in the Qqti 
system, although the latter is suppressed by extra powers of pion momentum. Before the 
three-pion threshold the amplitude is completely characterized by the pion-f2o scattering 
phases which are different in the | and § states, <5i/ 2 (a;) and 5 3 / 2 (w): 

' ((n n) 3 \A\n (0)) oc \A\e i5 ^ (97) 



2Mn 

for arbitrary operator A. These phases, strictly speaking, should be included into tl/ 2 (u;) 
and r 3/ / 2 (c<;). The FSI phases also depend on the isospin state, l= \ or I = |, but we now 
abstract from the light flavor symmetries. 

The scattering phases emerge via iterations of the f2o-pion interaction along with renor- 
malization of the coupling constant, and they must vanish in the limit of small pion four- 
momentum due to the pseudoscalar nature of the interaction. This fixes the physical value 
of the coupling as its value renormalized at the threshold. We identify this coupling with 



the one in the 'bare' chiral Lagrangian Eq. (92). Although the higher-order terms in pion 



momentum, including both scattering phases, can formally be obtained in perturbation the- 



ory in gs*B-K from the Yukawa-type Lagrangian in Eq. (92), this would make no sense for a 
number of evident physical reasons. What is relevant here is that even a naive perturbative 
expansion yields different scattering phases 5i/ 2 (w) and 63/2(00)- 

The net effect of FSI is to introduce an effective phase e l5i ^> in Tj and to replace 



the threshold gs*B-K in Eq. (96) by an energy-dependent coupling. In principle, the phases 
together with the couplings are constrained by analyticity of the amplitude and its unitarity 
property which are restrictive before the higher thresholds open. However, these may not fix 
the amplitude completely since the resulting relations are not local in energy and depend 
on the multiparticle domain. Likewise, the solution generically admits resonances in a 
particular channel whose number, position in energy and residues may vary. 

We shall neglect all such effects in what follows, assuming that the bulk of them is 
included in the resonance contributions, and subtracting the latter largely results in an 
effective cutoff of the soft-pion amplitudes at a certain scale near or below the lowest 
resonance. In particular, we attribute the difference between the \ and | channels to 
the resonant states. 

Next on our list is L = 1 producing the true ^ radial and | states; in our treatment 
they will again appear along with the | + '.D-wave' states. For the heavy quark expansion 
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we need to calculate the matrix elements for the operator with two derivatives, 



' (n 7r|QiDiiAQ|fto(0)), (98) 



2M Qo 

contracted with Sji, e m ji and (S m jS n i + S m i5 n j— \5 mn 5ji) to yield the spin-0, spin-1 and spin-2 
operators, respectively. 

To evaluate the pole diagrams for them we need the couplings analogous to Eq. ( 93 ) , at 
zero momentum transfer. These expectation values, however, are required over the heavy 
hadron moving with a small momentum of order u> ~ /ihadr- They are related by Lorentz 
invariance to the generic rest frame matrix elements with up to two spatial derivatives. 
Namely, the matrix elements of a product of any number of full covariant derivatives TriQV lx + 
7r M form a Lorentz tensor of the corresponding rank, and its value in an arbitrary frame is 
obtained, regardless of tuq, from the rest frame components by the corresponding Lorentz 
transformation. In the case at hand we get 

-^—(Q (p)\Qn 2 Q\Q (p)) = (p 2 +£)*l%, 

1 u 2 
——(QoipjlQiTCjTit-Ti^QlQoip)) = - ^l ajl %, 

- (S m jS n i + S ml 5 nj - pmndji) (Q (p) | QlTjlTiQl Q (p)) = (2p m p n - \b mn p 2 ) tfj^o- (99) 



2M, 



n 

Additional p-independent terms are absent since no spin-2 or higher current can be con- 
structed with a spin-| particle at rest. Similar relations hold for any higher order product 
as well; basically, the result for a non-zero momentum is obtained by incrementing each 7f 
by p: 

(Q (p)\Q7r jl ...7r jk Q\Q (p)) = (Q (p = 0)\Q(n J1 +p ] J...(n, k +pJQ\Q (p = 0)). (100) 

Otherwise the calculation of the transition amplitude proceeds exactly like for P-waves and 
yields 

-i-(O 7r|07r^g|O (0)) = g B * B7T - ^(^(k^-ha,) + (afykjk) (101) 

ZMq UJ V / 

The amplitude above calculated for spin-0 and spin-1 operators has a nontrivial 'diago- 
nal' Qq matrix element at rest, with the diagram a) not vanishing. The 'diagonal' piece of 
the matrix elements in Eqs. (99) (in this case it is the value at p = 0) is independent of the 
momentum. The propagators in diagrams a) and b) have opposite sign; this contribution 
then enters universally as a commutator with the pion interaction Hamiltonian, in accord 
with stationary perturbation theory: 

_ [6H,H] m n nn2 v 

Amn ~ {E m - En y ■ (102) 
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For instance, the kinetic expectation value always drops out, but a contribution remains 
proportional to the chromomagnetic interaction /j, G which does not commute with the spin- 
dependent pion vertex. This yields the most IR-singular contribution for a soft pion due to 



the double pole in Eq. (102). 



Strictly speaking, the pole diagrams with virtual f2o - they have a pole at u = (k is 
assumed to be fixed, and the overall k from the pion interaction is factored out ) - does not 
describe the pion emission amplitude completely. It contains a piece finite at u = coming 
from the contact interactions, see Fig. |T2fc . These contact vertices assume some values in 
QCD; we only know they are proportional to k since the heavy quark operators we consider 
are chirally sterile so that their (light flavor) axial charge vanishes: 



lim (Qon^k^QiDjiDiQino) 



1 



d*x (Q ir a \[QiDjiDiQ(0), J^ a \o,x)]\Qo) = 



(103) 

Otherwise they are arbitrary, and we can only parameterize them by four constants of 
dimension mass: 



1 



(n Tr\QiDjiDiQ\n ) = a ^(fc^-foo^o + p ie jZm A; m *5*Q + 7 5ji¥ iak^ + 



9B*Btt 



(104) 



plus higher terms in k. Such terms are induced, in particular, by intermediate excited heavy 



quark meson resonance states propagating in diagrams Fig. 12, a),b), with their relative size 
depending on spin. Unequal meson-to-pion couplings in Sect. 5.2 correspond to particular 
values of the constants a to 7; based on the QCD sum rule analysis Ref. [22] they appear 
large. While these terms are parametrically smaller than the pole amplitudes for small pion 
momenta, in general they contribute, especially the quantities regular in the chiral limit. 
There are also terms next order in Pn/Hh&dr- They would have the same scaling as the 



last term in Eq. (101). Nevertheless, we shall neglect them in what follows when address 



the pion loop effects proper: we relegate the corresponding contributions to the effects 
of resonances. This parallels the pol e-dominance assumption employed in the calculation 
of the D(*>ir amplitudes in Sect. 5.2 however the Z-diagrams included in the relativistic 
propagators of intermediate B* and D* mesons also induce such terms. If an argument can 
be put forward that one should retain only the D* or B* in the normal diagrams, there is 
no physical reason to exclude the higher states in the Z-diagrams: they all have the same 
large virtuality. Therefore, it may be natural to assume that in aggregate the Z-induced 
contact terms are suppressed or vanish in the nonrelativistic expansion. 
The result for the pion loop in R{uS) takes the following form: 
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+ -^-—^-{5ik&ji - 5ki<Jij + SijUki - SjiVik) + -^-{5ij5ki + S ik Sji + 5u5jk) 
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This can be decomposed into the invariant structures; introducing the common factor 
we get 



Pp (W) 



pf\u) 



W 
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8p 2 G \k\ 5 
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3w 2 

1 KPg) 
9uj 2 



2 ^k\ 3 



(106) 



The factorization properties stated in Sect. 6.1.1 are manifest here. 

For practical applications it is important to include in the amplitude the next-to-leading 



'contact' terms described by a, (3 and 7 in Eq. (104); they can actually be derived in a model- 



independent way and do not depend on the form of the resonance ansatz. Calculating the 
full RJj kl (u) with them is straightforward, but we pospone it to future work. 

The above calculation of the pion loop can be readily generalized to higher-dimensional 
heavy quark momentum operators, including the case of different number of derivatives in 
the two vertices. The expressions are particularly simple where no timelike momentum or 
antisymmetric spatial indices are involved; then the corresponding vertices simply amount 
to products of pion momentum and the resulting integrals can be simply calculated. 

We now return to our practical need. We are concerned with only the zero-momentum 
correlators of the l/m^-terms in the heavy quark Lagrangian, and we use the above derived 
spectral densities to estimate 



5PI 



and 



Sp 



~3 



0.02 GeV 3 



0.03 GeV 



-G 



0.065 GeV 3 , 8p% = 5p\ « 0.03 GeV 3 , 



8PI 



G 



0.1 GeV 



5p 3 s = 8p 6 A « 0.07 GeV" 



(107) 
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where we have assumed the upper cutoff at oj = 700 MeV and have adopted gs'B-K = 
4GeV _1 . The loops with charged and neutral pions are included, but not the kaon and 77 



contributions. Of special interest is the 'hyperfine' combination of Eq. (78) for which we 
get 

K-pIg-Pa)~ -0-09 GeV 3 ; (109) 

as expected, the pion contribution is suppressed, but it is negative and, taken at face value, 
would strengthen the lower bound on the resonant contributions. This is due to the positive 
sign of <5p 3 G , opposite to the BPS regime. 

With the same choice for g%* b-k we would have 



^0.015(^) 2 , 
^^(of^) 4 °- 006GeV2 



5A 



U.5GcV 

S Pd ~ (a 



5 GeV 



) 3 10 MeV, 

) 5 0.0025 GeV 3 



(110) 
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where we have anticipated a lower effective cutoff in the P-wave channel (the spin-triplet 
counterparts are not affected). We see that generally the pion loop contributes only a 
small fraction of the nonlocal correlators and a tiny amount of the local expectation values. 
A possible exception are quantities vanishing in the BPS limit where the pion loop may 
constitute a significant part of the deviation. 

The first calculation [33] of the chiral correction to the formfactor to order 1 / tuq ac- 
counted for the terms proportional to (pg) 2 from the HQS-breaking masses in the meson 
propagators and picked up only the p| and p\ pieces in Eq. (108), yet leaves out pl n and 
p^ G . We see, however that the latter contribution has the same size even though it is not 
singular in the chiral limit. 

One can combine the above description of the nonresonant D^ir states with the method 
of Sect. 6.3.2 for the inclusive yield, to get an estimate of the total continuum contribution 
to the inclusive yield. To that end one only needs to integrate the expressions in Eq. (106) 
over the mass of the Qoii state; due to the strong phase-space suppression the integral is 
effectively cut off at relatively soft pions. In this way we arrive at the yield in the ball park 
of 1% of T s \. This, however, refers to only the specified j p of the D^*\ states. The P-wave 
continuum is an independent channel. 



6.5 Nonfactorizable contributions to higher-dimensional local ex- 
pectation values 

The size of the expectation values of local heavy quark operators of D = 7 and D = 8 is 
important to estimate the impact of higher-order power corrections both in beauty and 
charm, to assess the accuracy of the OPE predictions and to study the convergence of 
the OPE series. The relevant operators here are those which emerge in the calculation of 
power corrections at tree level; they are sometimes called 'color-through' operators. Ref. [3J 
illustrated their effect in inclusive B decays for the semileptonic b — > c transitions and for 
B — >X s +j used in measuring \V c b\. They have also appeared in Sect. |4] in our analysis of 
B^rD*iv at zero recoil. 

In order to estimate the significance of the expectation values a ground-state factoriza- 
tion method has been devised in Ref. [3j which contains a derivation of the formalism and 
the explicit expressions for the factorization contributions to all nine dimension-7, m^g, 
and eighteen dimension-8, r!_ 18 , P-meson expectation values. Using the intermediate state 
saturation representation of Ref. [3J together with the relations elaborated in the previous 
sections we can now supplement the ground-state factorization values with the contribu- 
tions from the excited states: the analysis of hyperfine splitting allows to quantify their 
effect. This enables us to assess the accuracy of the factorization ansatz, the potential 
scale of the corrections to factorization and, ultimately, may elucidate the pattern of the 
higher-order effects in a more quantitative manner. 

Once again we start with infinitely heavy spinless quarks. The intermediate state rep- 
resentation for the operators with four spatial derivatives reads 

^-(^o|Q t 7Ti7r i 7rfc7ri(5(0)|Qo) = J dw R ijk i(uj), (111) 



45 



with Rijki(u) introduced in Eq. (47). The factorization contribution is located at u = 0; 
the upper limit of integration over u is determined by the normalization point assumed for 
the operator. The non-factorized pieces take a different form for each class of intermediate 
states. 

Contracting indices in Rij k[ (u) and multiplying the tensor, for the spin-triplet operators, 
by the spin matrix with an appropriate index, we obtain for the actual B mesons 

*„, = yp + yp + 2,<r> rt» = -i^p - P p 

^ mi = l P P + P p-ypP- P p-6 P ^ 

f 2 (i+) 1 fs+s 8 ri-h (2+-, 

3 _ 4 d m 7 — ~-^Ppg - ^Pf g 

5^m 8 = - 8p2 +) 

Tng—-—Ppg + Pg - —p f - Pfg -^Pg + ®P ' " { iiZ ) 

In the above equations the integration over u> is assumed, it has not been shown explicitly 
for compactness. We remind that m-i and are given by the fourth moment of the stan- 
dard SV (P-wave) structure functions and have no conventional ground-state factorizable 
contributions. 

For numeric estimates one simply considers the contributions of the individual multiplets 



of the excited states using their spectral densities in Eq. (|61l). The contribution of \ + is 



2 

obvious beforehand: it follows the factorizable one, see Ref. [3| , and only requires to replace 
(p 2 ) 2 by P 2 , (pg) 2 by G 2 and pip 2 G by PG. The effect of the higher-spin states has a 
different structure. Since the basis {mi — wig} has been selected arbitrarily, the impact of 
nonfactorizable contributions should be gauged in specific cases; this is easily done based 



on Eqs. (112). 



To quantify the overall scale of the effect we consider here three representative com- 
binations Mi,M2,Ms corresponding to the expectation values bif 2 if 2 b, b{aB){BB)b and 
—b(aB)n 2 b, respectively: 

Mi = mi + Jm 3 + ^m 4 , M 2 = -m 3 , M 3 = - Jm 8 , (113) 
2 6 o 



for which we have 



Mi = (/i 2 ) 2 + J dupf\u) 

M 2 = 2 -(p 2 G ) 2 + j dco (lpf\u)+pf\u)) 



M 3 = plp 2 + I dup { p p(u). (114) 



Numerically the corrections depend to some extent on the ratio of the | and | contri- 
butions and on P/G in the latter. Taking, for instance, e ra d ~ 700 MeV, J duo p g 2 (u) 
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6.3) and using the hyperfine constraint 



/ dujp p f\io) « / diop^ + \uj) (P*sG, see Sect. 
Eq. (80) with krj— 0.2 we obtain 

Mi w 0.206^ + 0.17 GeV;.^ 
M 2 « 0.08 GeVf ac t + 0.28 GeV^t 

M 3 « 0.15 GeVf act + 0.17GeVt fact . (115) 

A more definite value is obtained for the special combination 

-|m 3 - |m 6 - \m s 

which has the same structure as the hyperfine constraint; here we get 

- \m z - |m 6 - > 8 ~ -§(/4) 2 + 2/4^ - £ r ad(p^ G +pi) « (0.23 + 0.32 nf ) GeV 4 . (116) 

From this brief comparison we conclude that the factorization ansatz generally provides no 
more than a reasonable starting approximation for the expectation values not affected by 
cancellations. 

The hyperfine constraint and the approximations we complemented it with have nothing 
to say about the contribution of spin-| states which require a different theoretical input. 
Positivity and various relations between different nonfactorizable contributions are implicit 
in Eqs. (112); these constraints should, in principle, be applied only after the contributions 
from the non-resonant continuum are subtracted. 

The D = 8 operators with five derivatives were found to contribute at a lower level to 
the inclusive moments in B decays [3]; their precise expectation values are therefore less 
important in practice. Nevertheless, in this case it would also be useful to have at least a 
crude estimate of the potential error in the factorization ansatz, to be more confident in 
the assessment of the impact of 1/mn terms. 

The nonfactorizable effects for the D = 8 operators can be analyzed along the same lines 
as the D = 7 operators, identifying iDq adjacent to the intermediate state in question with 
—e. Most notably, the excited states contribute to a few combinations of t% which vanish 
in the ground-state factorization: 

r 5 = - U p { ) + \u) « -0.13 GeV 5 r 15 = -u- PF f ] (w) « -0.06 GeV 5 



2 

r 6 -n = - ui ( -pi 5 ] + -pV ] + ^py ] + -/a" ] + 3p ( * ! 



r 7 =~co ^-pi 2 } - -pi 2 ) + -p> } - -pi 2 } + 3p<* ) 

ri6 _ -UJ\ --Ppg - - Pg - — Pf - ~P fg + - P g + 3P 



, 2 1 3 (| + ), 1 (f + ), 1 , o (D 

r 17 -r 8 = - u> [ - Ppg - -p g - —p f + -p fg + - Pg + 3p^ 



r^-^lpP-^pp-lpf^Sp^). (117) 
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In the above equations the integration over cu is again assumed; the numeric estimates for 



rs and ri5 are obtained under the same assumptions as Eqs. (115). 



The most general analysis of the nonfactorizable corrections for the other operators 
(except for r\ and r 8 given by the fifth moment of the generalized SV structure functions) 
requires considering 



1 

2^ 



d 3 x / dxo e 



2M r 



-(ftolQWjQOr) Q^ k n mQ(0)\B) 



;ii8) 



which is an analogue of the tensor spectral density Rijki(co) in Eq. (47). One of the operators 



now includes an extra time derivative. The corresponding decomposition is lengthier than 



Eq. (48 ) since there is no symmetry between pairs of indices. The factorization of Eqs. (61 ) 



is modified for the new invariant structures: along with the transition matrix elements in 



Eqs. (58), (59) and (60) we need to introduce the similar ones P, G, f,g and h for the 



operators QttittqttiQ (it is assumed that ir acts on the right which matters for non-diagonal 



matrix elements). The analogue of Eq. (Ill ) for the operators with five deriva tives in t erms 
of Rijki{uj) holds and the general relations for the remaining similar to Eqs. (112) or (114) 
can readily be derived. 

The hyperfine splitting constraint cannot, however be directly applied to Rijki(w) and 
the corresponding tilded residues remain largely unconstrained even in the single excited 
multiplet approximation. The BPS approximation in this case yields P = G and f — g, but 
it is not too helpful. Therefore we do not quote here the corresponding expressions. 

It is nevertheless possible to get a rough estimate by making the assumption that the 
first P-wave excitation(s) approximately saturate, as an intermediate state, the transition 
amplitudes into the radial or the P-wave states: 



(119) 



3+ Qr 5- 
2 ul 2 



states under consideration. This can also be 



where p generically refers to the | 
regarded as an approximate relation obtained by truncating the complete representation 



7Tj|^ ) = ^ \f%e r 



M | v (m)\ . 



'3/2 



m n 

(120) 

after the lowest P-wave families. Such an assumption seems to work satisfactorily for the 
transition between the ground states, yet may be expected to degrade with higher initial 
and/or final states. Identifying e\/2 and 63/2 with e this would yield 



P ss -eP, G 
for individual residues, and 



-eG, 



9 



-eg, 



h 



-eh 



~( L+ ) ~ ( L+ ) 



for all invariant tensor structures with 1 = 1 and Z = 3, and, most generally, 

Rijki{u) w -e • Rijkiiui). 



(121) 
(122) 

(123) 
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Adopting for orientation such an approximation we obtain 



cP f r 2 ^ -ep { ) ] « -0.07 GeV 5 

* - i Qpf » + \?f > + ipf » + \ P f > + 3 P <! + ) 

5 nf r 9 ^5 nf r 10 « - £pJj +) « -0.07 GeV 5 

^nf „™f „ ~A (| + ) 3 (§+) 1 (| + ) (f + ) 

(3 rn«<5 r i2 I 6 ~ 40 ^ _ 8 



where the integration over u; is assumed similar to Eqs. (112). Taking e ~ 400 MeV the 
numerical estimates for the relevant combinations of the expectation values are then as 
straightforward as those for the D = 7 operators. The nonfactorizable corrections to the 
typical non-suppressed expectation values are of the order of 50 to 100%. 

Bearing in mind the dependence of all the contributions on a few poorly known hadronic 
parameters we cannot regard the calculation of general nonfactorizable effects accurate. 
They should be used primarily to assess the potential scale of the corrections to factorization 
and to clarify the expected sign pattern. The corrections associated with the | + states, for 
instance in {n 2 7r 2 ) or in (7? 2 7r 7? 2 ) are largely unconstrained since only physics related to the 
2 and | have been considered. 

The formalism employed in this subsection could be used together with the results of 



Sec. 6.4 to find the nonfactorizable contributions due to non-resonant states composed of the 
ground-state multiplet and a pion. However, we expect these contributions to be relatively 
small, as they come with higher powers of the excitation energy which is lower for the soft 
pion continuum than for the principal resonances, and therefore we do not consider them 
here. 



6.5.1 Ground-state factorization and N c 

The expectation values like (Qa-BxBQ) or (Qa-ExEQ) are possible due to the non-Abelian 
nature of QCD; in the factorization approximation they are proportional to |(/ic) 2 or to 
—e 2 fiQ. Such expectation values must vanish, on the other hand, in the bound states of 
Abelian theories like QED. In QCD proper one can consider similar time correlators for the 
usual magnetic and/or electric fields; they would enter, for instance, the electromagnetic 
corrections. The correlators can be decomposed into the same set of invariant spectral den- 
sities with different residues. For such Abelian fields certain expectation values must vanish 
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reflecting the commutativity of the different components of the Abelian field strength. For 
instance, this is the case for the counterpart of the p\ structure: 



dujp A 



Abel 



uj>0 



dui 



oj>0 



2 (| + ) 

3 "g, Abel 



f 3+i 
2^S,Abel 



2(fi 



G,Abel) 



(125) 



must hold for the Abelian analogies of the spectral densities and of the 5-meson local 
expectation values. This may be reminiscent of Eqs. (67), but is more general. Such a 



sum rule shows that the ground-state saturation itself may not be a universally applicable 



approximation. Eq. (125) may even be regarded as an indication of the importance of the 

3" 
2 



state. 



In actual QCD such local antisymmetric products of the guon field strength do not 
need to vanish; all the considered correlators and the factorized pieces count as constants 
in the large- N c limit. One may think that the ground-state factorization approximation is 
generally representative at not too small N c . 



6.5.2 Perturbative normalization point dependence 



Having at our disposal the perturbative heavy quark spectral functions of Eqs. (A. 22) 



and (A.17) we can easily incorporate the leading powerlike mixing for the operators to 
order a s and to any BLM order using the preceding analysis. The corrections to the 
expectation values are obtained from Eqs. (112), (117), etc.; the £>-meson spectral densities 



themselves are given by Eqs. (67). This gives the one-loop renormalization-scale evolution 



of the expectation values, except for the scale-dependence of the factorizable contributions 
themselves, like (nK/i)) 2 , which in practice may be significant. 

As a typical example, the one-loop and two-loop BLM piece in the combination M2 in 



Eqs. (114) is 



M 2 per V) 



r <*s(M) , 4 

Oj? jJL 



71 



1 + 



PqOCs 

2tt 



, M 29 

In 1 

2fi 12 



(126) 



with M denoting the normalization scale for a s (in the MS scheme). It is obtained using 



Eqs. (A. 25). The first-order correction to Mi vanishes and the second-order term is negative 
(although quite suppressed)]^] Of course, the perturbative calculation is meaningful only 
for not too low values of \i. At a numeric value of the strong coupling such estimates can be 
used to gauge the relative importance of the scale-dependence effects. The one-loop term 
with fixed a., = 0.3 would yield 

M 2 pert (0.7GeV) « 0.03 GeV 4 , 

and is relatively small for \x between 0.7 and 1 GeV. 

To extend this to the complete set of D = 8 operators including those in Eqs. (124) an 



additional class of the spectral functions would be needed where one of the pair products 



This shows a typical problem of applying naive non-Abelization to cases where the order-a s effect is 
absent; the full 0(a 2 s ) contribution to is, of course positive, paralleling the similar term in the Abelian 
theory. 
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71^71"; is replaced by Ttk^oTth °f- Eq. (118). The one-loop answer for it is simple, because the 



transition amplitude for such operators into a state with a single extra gluon amounts to 
— oj times the amplitude for the corresponding operator without extra n^. 

{Qg\QK k ^iQ\Q) = -ko(Qg\Qn k mQ\Q) + 0(g 3 s ). 



In other words, e = to in terms of Eqs. (124) to this accuracy. This exactly parallels the 
result for the P-wave SV operators. In this way both 0(a s ) and the higher BLM corrections 
to the mixing are readily obtained for the D = 8 operators alongside the power mixing for 
D = 7. 



7 Discussion 

7.1 On a resummation of the l/m k c corrections 

The accuracy in the estimate of .F(l) is limited, in particular, by significant higher-order 
power corrections in l/m c . We mention here a possibility to consider all these potentially 
dangerous corrections together and, therefore, in a certain sense, to resum them. The price 
to pay is the appearance of a limited number of new hadronic expectation values. 

The idea is to apply the OPE to the reversed zero-recoil transition D* — > B instead of 
B^-D*. At first glance this amounts only to exchanging c and b in Eq. and taking the 
expectation values over the vector rather than pseudoscalar statej^] However, the axial cb 
current produces not only B out of D*, but also B* . Therefore a scattering amplitude of 
two arbitrary spatial components j, I should instead be considered, and the corresponding 
indices contracted with the polarizations of D*, 

T zr (g ) = Jd 3 xJ dx e -^^-(D*\l iTc ljl5 b{x) b lkl5 c(0)\Dl) . (127) 

The analogue of Eq. ^ will still contain terms with powers of l/m c . However, they 
only come from the nonrelativistic expansion of the full-QCD charm quark operators cO^c, 
starting with the leading cc. Consequently, the full set of the l/m k c corrections originate 
from the heavy charm expansion of the finite-m c expectation value over actual D* states, 
(D*\cc\D*). The latter is a physical quantity and can in principle be measured on the 
lattice. Similarly, for a given power I of 1/mb all the terms l/m l b l/m k come from the the 
heavy charm expansion of the finite-m c expectation value of the corresponding cc operator 
with I derivatives. Likewise, the l/m c power effects in the inelastic transition amplitudes 
combine to yield directly the non-diagonal transition probabilities for the finite-m c charm 
states. The nontrivial explicit OPE corrections to the corresponding sum rule may only 
depend on powers of l/mt, since they come from the dynamic expansion of the intermediate- 
quark propagator. For what concerns the perturbative corrections, they are the same as in 
the sum rule Eq. ([I]), modulo the m c <H-m{, replacement. 

9 The fact that D* is not a stable particle is not an issue since the width of D* is extremely small. To 
completely bypass the complication one may simply assume that the pion mass is a few MeV larger than 
it is in reality; the properly defined B— >D* formfactor may not depend on this. 
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An attractive element of such an approach is that the m c -dependence of these matrix 
elements must be regular in the whole m c range and it should be possible to construct 
accurate interpolating functions, which would replace a large number of different coefficients 
appearing in the expansion in l/m c . Theory- wise, the regularity puts constraints on the 
higher-order behavior of the l/m c series and suggests that the series is sign-alternating, 
implying numerical cancellations between successive orders for actual charm mass. 

The series in l/m& is not resummed in this approach. As was discussed in Sect, [ij 
numerically these corrections can be discarded to a good approximation, and retaining the 
known leading terms for them would be sufficient for all practical purposes. 

Guided by the quantum-mechanical interpretation of the sum rules it is not difficult 
to verify that the leading, 1/wig corrections are identical to those in the direct approach, 
separately for the local OPE piece and for the inelastic contribution. Therefore the actual 
difference between the expansions would appear when the higher-order power corrections 
are addressed. 



7.2 Vector formfactor in B— >D transitions 

The present study focused on the B^-D* decay mediated by the axial current. A similar 
analysis may be applied to the B^-D transitions where only the vector current contributes. 
Two different aspects can be addressed here. 

The direct zero-recoil vector-current analogue Fd of J^n* is related to the matrix element 
(D\cyob\B) and does not determine the semileptonic decay rate near zero recoil for massless 
leptons, unlike B — > D*. It can be measured in the decay B — > D tv t whose amplitude is 
proportional to m T at q = 0, and this may represent an interesting opportunity for a new 
generation Super-B facility. 

The more conventional decays B — > D iv with nearly massless leptons are P-wave at 
small recoil and are more difficult to measure in this corner of the phase space. Of the two 
general vector formfactors f+(q 2 ) and f-(q 2 ) the latter does not contribute for massless lep- 
tons; the former, on the other hand, depends on both the time and the spatial components 
of the current (see, e.g. Ref. [20]). The spatial component assumes a change of the heavy 
meson velocity, is not related to a conserved Noether current in the heavy quark limit, and 
generally suffers from linear power corrections 0(l/rriQ). This fact fed, for a long time, a 
theoretical prejudice against the precision evaluation of the corresponding formfactor 

+ M B + M D J + Ky B D! 1 



It was nevertheless argued later [211], based on the BPS expansion, that in fact the power 
corrections in J r + are smaller, and may even enjoy a better numeric control than in Td*'- 

JF + = 1.04 ± 0.01 pcrt ± 0.01 power . (128) 

The analysis developed in the present paper can be applied to both formfactors, Fd and 
J-+. For Fd the required modifications are minimal. The case of J-+ is somewhat different 
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both technically and conceptually; in particular, the physical interpretation is different, and 
there is no simple probabilistic interpretation that would mean positivity already for the 
1/mg corrections. The positivity holds for the leading l/mq power correction, however it 
can simply be regarded as known within the required numeric precision. 

The analysis of the perturbative corrections follows closely that of Sect. [3] in the case 
of Fjj, with corrections generally smaller. This kind of analysis cannot be directly applied 
to due to subtleties at nonvanishing recoil; the treatment is more complicated here and 
has to be analyzed ad hoc [20J. 

For what concerns the power corrections, we do not expect improvements in either 
cases. The main reason is that the power corrections (in the kinetic scheme we consistently 
use) are numerically small to start with, since they all vanish in the exact BPS limit. 
Moreover, to any order in 1/mg the terms are of the second order in the deviation from the 
BPS limit [20] , an analogue of the Ademollo-Gatto theorem [34J which applies to the BPS 
expansion for both Fp and T+. As a consequence, any numeric result strongly depends 
on the degree of proximity of the actual-QCD dynamics in B mesons to the BPS limit, for 
instance, on the excess of fi^ over \j? G . This appears as a significant cancellation between 
terms belonging to different spin structures. It is therefore difficult to expect an increase 
in defendable accuracy in Fd and J-+ unless the numerical aspects of the BPS breaking are 
experimentally scrutinized. Once this aspect of the strong dynamics is studied, we would 
have more theory constraints to improve the accuracy of the nonperturbative predictions, 
in particular if the BPS regime turns out a good starting approximation. 



7.3 Lattice determination of Tr>* 

The recent PDG policy has been to rely solely on the lattice evaluation of J-"d* for the 
exclusive extraction of \V c b\ from the B — > D*lv differential rate extrapolated to the no- 
recoil kinematics. The lattice values for have always been on the higher side, well 
above 0.9 and carried small error bars, in particular since unquenched simulations were first 
employed [35] : 

T D *(\) = 0.924 ± 0.012 ± 0.019. (129) 

An update of this result was presented by the FNAL-MILC collaboration |36j after our first 
publication pQ. It has a lower central value rj 

T D *{1) = 0.902 ±0.005 ±0.016 (130) 

and is closer to our number. The other recent lattice result is based on a quenched simulation 
by the Tor Vergata group [37], •7 7 d*(1) = 0.924 ± 0.008 ± 0.005 and has an even smaller 
nominal error. 

Confronting our evaluation of J r o*(l) with the lattice ones, we should first emphasize 
that the latter are not direct calculations of this formfactor. The lattice theory with heavy 
quarks and continuum QCD are two different theories, and there is no limit at mq a ~ 1 (a 



10 



The last FNAL paper included the electroweak decay enhancement factor 1.007 into F; we have removed 



it in Eq. (130). We thank A. Kronfeld for the communication. 
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is the lattice spacing) where they would coincide nonperturbatively. Present simulations do 
not reach values of m c a below about 0.3-0.4. The fact that the two theories share the same 
heavy quark symmetry was emphasized as the key point behind the approach pursued by 
the FNAL group [38] • However, it is the power-suppressed deviations from the symmetry 
that matter in this case, and in principle they are different. 

This becomes transparent if one takes a closer look at the l/m,Q corrections. In QCD, 
as a consequence of Lorentz symmetry, the mass entering the nonrelativistic kinetic energy 
p 2 /2rriQ is the same rest-energy mass itlq controlling also the power-suppressed terms in 
the currents (from equations of motion and from the Foldy-Wouthuysen transformation), 
the role of the correlators of the subleading operators, etc. No symmetry, however, enforces 
their equality on the lattice, and all these terms are driven by different effective masses. 

The FNAL lattice approach appreciates this complication [39]. To handle it the heavy 
quark sector is modified by adding ad hoc power-suppressed terms allowing to change the 
corresponding effective masses; in turn this changes the power corrections. The outcome 
for the effects in .F(l) is then determined by the ad hoc constants which have to 

be specified through a matching. However, such a matching has only been performed at 
tree level. Furthermore, the matching was performed only to the leading power effects, 
corresponding to the l/m^ corrections in .F(l). All l/rriq and higher effects are therefore 
not under control Fj although it can be argued that these discretization effects are somewhat 
suppressed for charm, at low a. 

In view of these practical limitations the FNAL approach cannot be regarded a first- 
principle evaluation of the zero-recoil B — > D* formfactor in QCD. Some of the related 
potential biases are included in the error budget detailed in the publications. However, the 



error assignment may not be realistic. In fact, the recent value in Eq. (130) has a reduced 



discrepancy with the estimate Eq. (133); a more conservative treatment of the systematic 
errors would make it compatible with the central value of the present analysis. 

We also note that, comparing the earlier and the more recent FNAL lattice simulations, 
the group did not find a noticeable effect of the light quark unquenching. This differs from 
the estimated size of the chiral loop contributions discussed in Sect. |5.2| related to the non- 



resonant states with light dynamic pions, although there may be no formal contradiction. 

The above reservations apply to the lattice determination of the vector B-^D formfactor 
T+ as well. In this case larger corrections to the symmetry limit were found 

7+ = 1.074 ±0.018 ±0.015 



compared to the theoretical prediction in Eq. (|128|). Here, however, the disagreement is 
less significant than for Td* 



12 



11 Earlier FNAL evaluations [JD] claimed to extract the principal l/rriq terms, however they were deter- 
mined in an effective theory essentially different from the actual QCD. The later analyses considered only 
1/rrig corrections. 

12 More recently, the Tor Vergata group reported the value F + = 1.026 ± 0.017, based on quenched 
simulations [H] and a new preliminary value T+ = 1.058 ± 0.009 sta t [33] has been presented, where the 
systematic uncertainty still needs to be evaluated. 
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8 Conclusions 



The present study has been motivated by the need for an updated evaluation of the phe- 
nomenologically important B -4 D* semileptonic transition formfactor near the zero-recoil 
point, ^(1), that could account for the latest progress in heavy quark theory. The main 
numeric outcome has been reported in [1] and the details have been given here. 
Numerically we conclude that the unitarity upper bound for the formfactor is 

J"(l)<0.92 (131) 

assuming only positivity of the inelastic contributions. Including the soft D^ir continuum 
the bound becomes 

J r (l)< 0.90. (132) 

These numbers refer to values of /i^. close to its lower bound. The bounds on ^(1) become 
stronger if the actual ^ value is larger, which is more natural on theory grounds. Our 
analysis of the inelastic transitions into the low-lying channels incorporating the constraints 
following from the observed amount of the hyperfme splitting in B and D mesons allows us 
to go beyond the unitarity upper bound and to obtain an estimate. The actual value for 
J-"(l) comes out about 

J 7 (I) « 0.86 (133) 
at low values of it somewhat decreases at larger and/or p 3 D . 



The quoted numbers for the unitarity bounds Eqs. (131), (132) carry a theoretical 
uncertainty of about 1%. The estimated central value has a larger uncertainty, of around 
2%. We quote here the value literally obtained in our estimate (for low /i^, p 3 D ); it is not 
implied that this value peaks the expectation probability. 

Thus, J-"(l) in excess of 0.9 would be consistent with unitarity and the short-distance 
expansion of the QCD amplitude only under contrived assumptions. Values larger than 0.92 
should be viewed in violation of unitarity assuming that the conventional short-distance 
expansion in QCD works in the case of the zero-recoil scattering amplitude off heavy quarks. 

In earlier analyses of the power corrections to ^(l) the "wavefunction overlap" effect 
used to be uncertain and, essentially was only parameterized; in the language of the heavy 
quark sum rules it was 

/inel = X ■ A (134) 

with A the power corrections in the sum rule, cf. Eq. ([6]), setting the scale of the nonper- 
turbative effects in J-"(l). It was simply guessed following Refs. [10112] that x is somewhere 
between and 1 leading to an assumption x = 0.5±0.5. On the other hand we have linked 
/inei to measured hadronic parameters in the heavy mesons and found \ to be large, 

X ~ 1.3 4- 1.7. 

This is the main factor driving the prediction for J-"(l) down compared to earlier estimates, 
along with a shift due to the higher-order power corrections. 
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Since our conclusion appears in some conflict with the lattice results for .F(l), one may 
examine how robust this conclusion is. We found that the phenomenology of the heavy 
mesons suggests that the inelastic contributions, as well as the l/m^ and higher-order 
power corrections are numerically significant and that they lower the expected value for 
•P(l). Our derivations, of course, use the power expansion for charm, and this may be a 
vulnerable point for a precision prediction. 

We emphasize, however, that a scenario with smaller corrections from the nonlocal 
correlators or high orders would not be consistent. If the neglected higher-order effects 
in charm are relatively small, our analysis of the hyperfine splitting is reliable. If, on the 
contrary, the higher-order corrections in charm are too significant to affect the credibility 
of the analysis, this would imply a higher overall mass scale of nonperturbative QCD in 
heavy quarks. Then there would be no reason to expect small corrections to ^(1) either. 

Our numeric predictions involved a number of theoretical improvements. The first to 
mention concerned the calculation of the Wilsonian perturbative renormalization factor 
with a hard cutoff. We derived a general ansatz applicable to the one-loop level as well 
as to arbitrary BLM order, which yields the perturbative correction in the kinetic scheme 
with a full dependence on fi/rriQ. It appears to be different from the naive prescription for 
having a cut on the gluon momentum. For the zero-recoil transitions the difference emerges 
starting with the terms 0(1/i71q). Numerically it is important whenever the hard scale is 
determined by the charm mass. 

The other direction is the treatment of the inelastic contributions, along with a detailed 
analysis of the continuum soft-pion states in the context of the heavy quark expansion. 

We have presented a novel model-independent analysis of the transitions into the radially 
excited (or P-wave) states near the rest kinematics. So far the excited states considered 
were mainly the P-wave states. The important new phenomenological constraint comes 
from the hyperfine splitting in B and P: qualitatively, the latter tells us that the D = 3 
zero-momentum nonlocal correlators are numerically large in actual QCD, and this enhances 
the predicted size of the inelastic probabilities (the 'overlap deficit' in the formfactor) over 
the naive expectations. The analysis of the spin-averaged P and D meson mass difference 
supports the same conclusion, albeit with larger uncertainty. Using the hyperfine splitting, 
one can then set a lower bound on the inelastic contribution. The bound is very close to 
the value ii ne i assumes in the BPS limit where a single combination of the four general 
correlators determines both 7i ne i and the I/tuq dependence of the hyperfine splitting. 

A related implication of the hyperfine splitting analysis is the enhancement of the tran- 
sition amplitudes into the excited 'radial' charm states which must be dominated by the 
l/m c -suppressed terms rather than by the velocity-dependent component. This leads to an 
increased yield of the wide charm hadronic structures from the decays of the 'radial' states, 
and would eliminate the '§ > §' puzzle if the observed wide yields are dominated by them 
rather than by the | P-waves predicted to be suppressed. We emphasize that the P-wave 
states must be produced along with true radially excited mesons, and may even dominate. 
Some possibilities to distinguish them in experiment were discussed. 

Therefore, we have identified a link among three apparently unrelated physics points: 
the size of the hyperfine splitting in charm and beauty, the reduction in J-{1) through 
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the enhanced nonlocal power corrections, and the resolution of the '| > |' puzzle in the 
semileptonic B decays. 

On the theoretical side, yet another implication is manifest: we find significant cor- 
rections to the factorization for higher-dimension heavy-quark expectation values. Our 
model-independent analysis provides the basis for a dedicated account of the nonfactor- 
izable effects from higher orders in I/tuq [3J, in particular for the semileptonic -B-decay 
fits. 

With respect to the D^n states, we have expanded the treatment in a few aspects. Our 
approach allowed us to account for their effect in the analysis of J-"(l) beyond the leading 
order 1/rriQ, relying instead on the soft-pion approximation. The sub leading terms turned 
out significant. In particular, the heavy quark symmetry-breaking corrections in the heavy 
meson-to-pion couplings seem to yield the dominant effect. 

Theoretically, we have presented a consistent treatment of the D^tt states in the heavy 
quark approximation within the soft-pion approach. The decomposition of the P-wave into 
the |- and the |-components has been addressed before yet remained largely unknown. 
We briefly recapitulated it in Sect. 6.4 and extended it to radial/D-wave states, with the 
decomposition into the corresponding ^ -, | - and | -channels. This made it explicit that 
earlier studies of the soft-pion corrections were incomplete: they accounted only for the 
most singular effect at the softest pion momentum, however not dominant in the typical 
configuration with ~/ihadr- Typically, the pion loops yield negligible contributions, with 
the notable exception of the contribution to the axial sum rule for ^(l). 

We conclude with the following remark. The pattern of the hyperfine splitting in B and 
D mesons, in particular its precise mass dependence, is important for a few different phe- 
nomena and draws novel qualitative conclusions for our understanding of the heavy meson 
states. The precision interpretation of the splitting, on the other hand, may potentially be 
hindered by higher-order effects in charm. A sufficiently accurate measurement at a dif- 
ferent heavy quark mass would radically improve the credibility of the hyperfine analysis. 
We have pointed out that a first principle lattice determination of the hyperfine splitting 
at better than 5% accuracy, if possible, can provide this information, and have discussed 
how it can be used for mesons either heavier or lighter than charm. 
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9 Appendices 



A BLM corrections and summation 

The technique that allows the computation and summation of BLM corrections of order 
/3qQ;™ +1 has been concisely reviewed in Ref. [8]; Ref. [7] focussed upon its application to 
Wilsonian OPE calculations. To BLM-dress a one-loop result one needs to evaluate the 
generic one-loop correction A\ 

A 1 • n M, • ... 

7T 

with a fictitious gluon mass A, Ai(X 2 ), so that the conventional A\ is Ai(0). The BLM 
series in terms of the MS coupling a s normalized at the arbitrary scale M takes the form 



.4 



BLM 



k=0 



Mo) 



a s {M) 
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2\k 



c 



2k + l 
n+1 
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ft) 



4-7T 



n+2 



X 



n=0 

dA 2 
A 2 



' M 2 5 
'"^ + 3 



-i n-2k 



A x (0) 



M 2 

M 2 +e -5/3 A 2 



-Ai(\ 2 



(A.l) 



with C denoting the binomial coefficients. Let us remind that the integral over A 2 here 
typically has two (or more) domains. The last term depending on A\ (A 2 ) is integrated from 
to the threshold value of the gluon mass if there is a threshold (in particular, it would be 
set by \x in the Wilsonian calculations), whereas the first term in the same brackets should 
always be integrated over all values of A 2 regardless of a cutoff or of kinematic details. 

When the infrared part is removed from the one-loop diagram to leave a genuinely short- 
distance correction, the corresponding Ai(A 2 ) is a real analytic function in the vicinity of 
zero. This allows to write an integral representation for the resummed series, 



A 



BLM 



1 + ^(0) 



a s {M) 



TV 



+ 



dt 



-co (l + ^(t-§)) 



4 

A) 



M 2 



M 2 ~A 2 QCD 



A 1 (0) - A 1 (-A 2 V ) 



A 1 (0)-\ ir -A 1 (e t M 2 : 
1 + e*"s 



(A.2) 



without any ambiguity associated with the last term; here 



A 2 QC D = M 2 e /WW). 



A 2 



QCD 



M e W") 3 



(A.3) 



B One-loop perturbative calculation with a Wilsonian cutoff 

In this Appendix we discuss aspects of the one-loop calculation of the leading Wilson 
coefficient £4. The reasoning is quite general and is applicable to other observables as well. 



58 



As in Sect. |3j we need to distinguish between em and Wilsonian fi, therefore we will deal 
explicitly with ^ ert (£M, A 4 ) as it appears in Eq. The dependence of ^^(em, A 4 ) on em is 
given by Eq. (|8]). Since the right-hand side of Eq. ^ must be in dependent, and £^ ert must 
satisfy the 'boundary condition' Eq. one way to determine ^^(^a/, A 4 ) f° r arbitrary 
values of /i and (in principle) any order of perturbation theory is to use the /i-dependence 
of the matrix elements and masses appearing in Eq. (J6]). 

Indeed, in perturbation theory the expectation values (B\Ok\B)^ in Eq. ^ typically 
depend on /z in a powerlike way, 

(a.4) 

d/i 

where dk is the dimension of Ok] the same applies to the heavy quark masses on which 
the Wilson coefficients generally depend. In this way the ^-dependence of ^^(emi A 4 ) 
is calculated explicitly while its /z-dependence emerges as an expansion in [ijvriQ which is 
necessarily truncated after a few terms. This is not a serious limitation if the hard OPE 
scale is 0(rrib), but in our case the charm mass sets the lower hard scale in the problem 
and the expansion in \ijmQ shows poor convergence. 

To overcome this drawback we have devised a method allowing to directly compute 
the leading Wilson coefficient £^ ert (£kf,/i) as a function of /i/mg, without an expansion, 
in the one-loop approximation. The method is readily generalized to higher-order BLM 
corrections. We have described its main points in Sect. [3j below we provide additional 
explanations. 

The idea behind the approach is that in one-loop calculations there is a simple connection 
between the normalization point of the heavy quark operators in the kinetic scheme and the 
hard cutoff on the gluon momentum k in the diagram. In order to preserve the analiticity 
and the unitarity of the Feynman integrals no limit on integrations over ko of the gluon is 
imposed in the kinetic scheme: the separation of scales is performed based on \k\. In this 
way one computes the one-loop £^ ert introducing an infrared cutoff on \k\, instead of taking 
the full integral d 4 /c. In the 6-quark static limit the step-function cutoff factor 

9{\k\-ti) (A.5) 

in the Feynman integrand yields precisely the normalization at the scale \x. The remaining 
part, the integral with 9(fi— \k\) constitutes the power-suppressed terms described by higher- 
dimension matrix elements in the OPE. In order to go beyond the static approximation 



certain modifications of the cutoff in Eq. (A.5) are required. 

Let us remind why £4 is related to rf A calculated with a cutoff. The reasoning is based 
on considering the OPE relations in an ensemble of gluons with the spatial momentum 
limited by /x in the h rest frame; the non-Abelian nature does not play a role at one-loop 
level. All these gluons can be considered soft, and they satisfy the OPE sum rule where 
the coefficients assume the tree level values. Having in mind how the sum rules are derived 
(for pedagogical reviews see Refs. [TTj [21]), the integration in the sum rule must run over 
all excitation energies, from to 00. However, in the soft gluon ensemble with \k\ < \i the 
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spectral density vanishes above the excitation energy 

£m{v) = A*+ \Jm 2 c +n 2 -m c ] (A.6) 

below £m(ji) it is the usual one-loop spectral density of QCD. At the same time, the operator 
expectation values in such an ensemble are just the one- loop QCD expectation values in 
the kinetic scheme normalized at fi. 

We now turn to the subtleties beyond the static limit, and focus on the soft contribution 
8rfj^ to be subtracted from the one-loop tja- 

Va^Va- ^ oft = VA - C F g 2 s J -0^. 6(n-\Z\) . . . (A.7) 

Here the ellipses denote the same propagators and vertices encountered in the calculation 
of t]A itself; k is the gluon momentum in the diagram. There are three one-loop diagrams 
- the vertex correction and the wavefunction renormalization for both b and c quark. 

Let us consider the vertex diagram of Fig. [2^, as an example. The integrand has a 
general structure 

d'7r numerator 



(2tt) 4 z (k%-k 2 +iO){m 2 b -{m b -k ) 2 +k 2 -iO){m 2 c -{m c -k ) 2 +k 2 -iO) ' 

At given k the integral over ko is convergent and is saturated at \ko\ ~ \k\; the tails at 
large \k \>rriQ contribute a power-suppressed piece. Since no cut on k Q is allowed, the 
integration over ko can be performed by closing the integration contour in the lower half- 
plane, see Fig. [2b. There are three pairs of poles in the ko plane, 



ko = ±\k\, ko = mb±y ml + k 2 , ko = m c ± y m 2 + k 2 . (A.9) 

With the standard Feynman prescription, blue contour a, the integration over ko results in 
the sum of the three residues corresponding to the above poles. For the first pole we have 
k < fM] however for the two other, distant, poles we have k >mQ 3> fi, regardless of the 
cutoff. Since the OPE generally corresponds to an expansion in all components of the gluon 
four-momentum, it is clear that the contributions to the integral associated with the distant 
(black) poles may not correctly describe the OPE power-suppressed terms. We recall that 
our goal is just to subtract the piece of the one-gluon loop correction to t\a associated with 
the terms which have already been included in the power-suppressed OPE. 

Indeed, it turns out that the OPE series for the soft piece correspond to the residue of 
only the 'near' pole at fen = |fc|, while the two other resides should be discarded. This means 
changing the bypass prescription for the two distant poles, —iO-^+iO, which moves the ko 
integration contour as shown by the green dashed line in Fig. [2b. 

The 1 / rriQ expansion in the Feynman diagrams leading to the OPE series is essentially 
the Taylor expansion of the heavy quark propagators in the integrand for small gluon four- 



momentum k. Eq. (A. 8) then becomes 



d 4 & numerator (k 2 ) n (k 2 ) m ^ 



(2n)H (]feg_fc2+i0) ^ (2m b k -i0) n+1 ^ (2m c k -t0y 
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the poles at ko = + zO are descendants of the nonrelativistic poles located on the left of 
ko = 0. To reproduce the nonrelativistic expansion we calculate the integral for each term 
closing the contour in the lower ko half-plane and picking up only the gluon pole ko = \ k\. 



This cannot be done for the original integral in Eq. ( A. 8 ): its value is not equal to the sum of 
the series, because large k beyond the convergence radius of expansion (A. 10) contribute. 
However, with the modified bypass prescription for the distant poles the integration contour 
in Fig. [2|d can be shrunk to a contour 7 (maroon) on which \ko \ <2m c holds everywhere. 
For the integral over contour 7 the series converges absolutely and uniformly, still embracing 
only the ko — \k\ pole. This proves that the sum of the series obtained integrating term by 



term Eq. (A. 10) gives the integral over fc of the original expression in Eq. (A. 8) yet with 
the modified bypass for the distant poles along contour (3. 

Let us make the following general remark. At first glance, the terms in the expansion 
of the Feynman diagrams and in the OPE series assume a somewhat different form: the 
expectation values of the operators in the kinetic scheme are given by three-dimensional 
integrals over spatial momentum of an on-shell gluon, k Q = \k\, being defined through the 
heavy quark structure functions: 

r\ 3 k 

-V(k) (A.ll) 



(2vr) 3 2|fc 

with V(k) a polynomial. It is clear, however, that the two representations have the same 
form once the k integration is performed closing the contour in the lower half-plane. That 
is why integration over k plays an important role in our reasoning. 

At this point it becomes transparent why the OPE series yield the expansion of the 
sole contribution of the near pole with the on-shell gluon: the heavy quark propagators 
appearing in the definition of the heavy quark expectation values in the effective theory 
are nonrelativistic (static) propagators which have a single pole. The second, distant, pole 
peculiar to relativistic particles is absent from them. Therefore, at any finite order in the 
OPE power expansion there is no contributions associated with the distant poles, k ^2mQ 



in Eq. (A.9). In this sense the distant singularities are related to the divergence of the 
expansion in kjrriQ rather than to the discontinuity of the individual terms. 

So far we have considered the vertex diagram. The other two Feynman diagrams with 
the wavef unction renormalization for external quark legs have the same structure; they only 
depend on a single quark mass mQ = m^ or m c , and the two pairs of the fermion propagator 
poles are degenerate. Consequently all the above reasoning applies to them as well. 

Let us note that taking \i = em the difference between fjf and Em becomes power- 



suppressed, the last term of Eq. (10) becomes of order 1/mg and can be neglected to the 
leading order fi 2 /rriQ. Therefore, it accounts for the recoil correction in the relation between 
\i and Em and becomes relevant where the terms 0(a s /i 3 /ttiq) are included. Its form must 
already be clear from the preceding derivation: in the soft gluon ensemble the emission of 
a gluon with energy uj yields an excitation energy e = us + \J m 2 -\-uj 2 — m c in the final state, 
and the spectral density is 



e) =w pcvt (E)9(fi'-£), fi' = fi+y/m 2 c + fi 2 -m c . (A.12) 
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This explains the integration limit in the last term in Eq.(lO) 



The above discussion of the one-loop corrections is directly extended to incorporate any 
higher-order BLM corrections as well, or even to perform the complete BLM-summation. 
As detailed in Appendix |X| the same analysis must be repeated for the diagrams with an 
arbitrary gluon mass A. This has been stated already in Sect. [3] where the related technical 
modifications were listed. The necessary one-loop expressions at non-zero A 2 are given in 
Appendix O 



C Details of the perturbative calculation 

In this Appendix we provide details omitted from the main text in Sect. [3] In the case of 
a massive gluon the one-loop inelastic perturbative spectral density determining the recoil 
correction is given by 



w peit (e) 
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F vr 12 M 3 (M — m c ) 2 (m b M 2 — MX 2 + m b (\ 2 — m 2 )) 2 
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+M 2 (4m 2 A 2 - Am b m 3 + 6A 



6m 2 b (m 2 c - 2Y 



(A.13) 



Here M = m c + e is the invariant mass of the final state. The one-loop zero-recoil axial 
current renormalization without a cutoff t\a = I + C-F^^a^A 2 ) + 0(a 2 ) is given by 



9y + 2yz + 24 + 9yz 2 9y 2 - 7y 2 z - Qy - 6yz + 18 + 18z 
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1 - A/yz 2 



A8y(l-z)^l-A/yz 2 



Here we have used y = A 2 /m 2 and z = m c /m b . This expression coincides with the one 
in Eq. (B3) of Ref. |Hj, where it was derived by a dispersion integral starting from the 
Euclidean calculation of Ref. [13]. 

The check of the //-independence of the OPE sum rule, Eq. (J6J), with the perturbative 
factor calculated according to Eq. (10) can be accomplished to an arbitrary BLM order 



at once: it suffices to establish it at a given value of the gluon mass. We demonstrate it 
here assuming e M = [i. The right-hand side of the sum rule depends on // through £^ ert (/x) 
and through the heavy quark expectation values. Moreover, since all the //-dependent com- 
ponents of the perturbative calculation, including the expectation values and the inelastic 
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integral, can be written as three-dimension integrals 



d 3 A; 



(2vr) 3 2£; 



00 - feo) • 



k = \^k 2 + X 2 



(A.15) 



with the same cutoff provided Em = P, we may simply check the cancellation at the integrand 
level at a given value of k 2 and A 2 . 

As the //-dependence of the individual contributions starts at 0{l/m 2 q), the first check 
is provided by the terms a s p 2 /rn^ [35]. At this order the power-suppressed component of 
the OPE part of the sum rule (|6j includes only the expectation values /i 2 and p? G which 
should be considered in the static limit; the latter then vanishes. The required expressions 
have been given in Refs. [8j [7] . In units of 



Cf9%s I W^k~ 



we get 



2 + 



A* 

k 2 ' 



1 



m 2 ' 3m c mj ' m 2 , J 2k\ \m c 

for ii 2 , and w pert , respectively. 

Finally, the subtracted soft piece of tja depends on p. The expression for Srj s ^ it (p) given 
in Eq. (18) must be expanded in k and A 2 : 
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(A.18) 



where the first square bracket upon integration gives the leading fi 2 /mq terms and the 
second line yields fi 3 /mg. Combining the above coefficients the /^-independence of the sum 
rule is verified at order fj^/mg. 

A far deeper check is encountered at the level of /i 3 /wig terms. At this order the proper 
prescription to calculate the soft virtual correction 8rf^ rt discussed in Appendix B becomes 
essential. At this level one also needs to include the higher-dimension Darwin expectation 
value as well as the l/m& effects in the kinetic and chromomagnetic expectation values. The 
latter are expressed in terms of the local p 3 D and p\ s and of the nonlocal expectation values. 
Perturbatively p\ s vanishes as do p^ G and p\. The subleading term in the continuum 



spectral density is also required and can be obtained directly expanding Eq. (A. 13). 



By virtue of the SV sum rules the perturbative component of the Darwin expectation 
value amounts to an integral of the same form (A.15), with the integrand given by the 



product of the integrand for /i 2 and ko, the excitation energy for a static quark: 



fcb(2+g: 



(A.19) 
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(b) 



Q 



Q 



Q 



Figure 13: One-gluon amplitude diagrams involved in the calculation of the 0(a s ) perturbative 
spectral densities in the static limit. The square block denotes the operator in question. 



Since p? G was defined in Eq. (20) to include only the magnetic field piece —o ■ B but not 
the lower-component term in the complete chromomagnetic operator, to order fi 3 /m,Q it 
acquires a contribution only from the nonlocal correlator p|, 



The value of p 2 is perturbatively corrected, instead, by the correlator p: 



,3 . 
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(A.20) 



(A.21) 



where overal integration (A. 16) is understood. 



To find one-loop 5 aa p^ n s or, more generally, all the perturbative spectral densities in 



Eqs. (66), (67) one needs to square the sum of the diagrams shown in Fig. 13. The answer 



is obtained immediately, in particular if a little trick [IB] is used, which eliminates two of 
the three diagrams, b and c, for the transition amplitude: 



7T 



2^1 



C 



3tt 



7T 



(A.22) 



The trick uses gauge invariance to say that one can use the simple gluon propagator 



n - 14 



(A.23) 



A 2 - k 2 

requiring to calculate only the spatial vertices; they come from the vertex emission but not 



from the heavy quark lines. Using Eqs. (66), (67) and (51) we obtain at once, for instance, 



2/„2 



3 

7V7T 



X 2 k 



SasPs 



(A.24) 



in the same units of Eq. ( |A~T6| . 

Finally, one needs the (9(p 3 /mg) term in the expansion of Srffi ; it is given in the second 
line of Eq. (A. 18). Expanding Eq. (A. 13) through the next-to-leading order, on one hand, 



and collecting all relevant terms from Eqs. (A. 18), (A. 19) and (A.24) with the explicit 



coefficients appearing in the sum rule, on the other hand, we arrive at the same integrand 
in the both sides of the sum rule (p| at 0(l/mg), for arbitrary value of A 2 . 
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For bookkeeping purposes we quote here the perturbatively calculated moments of the 
spectral densities po(u>), pi{u) and p^{oS), as well as those of the P-wave spectral density 
p p (o;), defined in analogy to Eq. (67) by the perturbative relation 



9 



see Eq.(44). Their expressions to order a s and Poaj. are: 

0(1 



u k doj po(co) 
uj h du pi(u) 

I 

oo k du p P (uj) 



Po /"A 2 P 4+k 



7T 



4+fc 
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2C F 



aJM) P 

7T 4 + k 
a s (M) p i+k 

7T 4 + k 
a s (M) p 2+k 



7T 



2 + k 



A + k 
1 + 

1 + 

1 + 



2 7T 

A) «£ 

2 7T 
2 7T 



ln M 13 1 

11 2^ + ~6~ + 4 + fc 



, M 71 

In 1 

2p 30 

, M 5 

In h - 

2p 3 



1 



A + k 
1 



(A.25) 



The moments of p F (uj) determine, beyond the IW slope, A, the kinetic and the Darwin 
expectation values, the higher expectation values ui2 and r\. 

For completeness, we give here also the perturbative dependence on p of p\ s associated 
with power mixing: 



^(-pls) = C A ^pl + C A ^ 1 -pl s + 0(al) 



(A.26) 



(the anomalous dimension of the spin-orbit operator coincides with that of the full chro- 
momagnetic one). This implies the relation between the extrapolated 'pole-scheme' value 
—p\ s an d the Wilsonian —Pisi^)'- 



~pls = -pls(p)-C A ^PP 2 G + 0(a 2 s ). 



(A.27) 



D Transitions to the 'radial' states and the inclusive yield 



The inelastic zero- recoil spectral density for the vector cjob current paralleling Eq. ( 52 ) is 

2 _(§ + ), 



2ni 



discTf)( 



i 



i 



2m c 2mi 



(i+) 

Pp ' (u)-2p p 2 g (U)+Pg 2 (W) 



(A.28) 



no transitions into | + occur. It is manifestly BPS-suppressed to the second order. 



As stated in Section 



6.3.2 



we evaluate the rate for a decay of the flo spin-^ heavy state 
into the corresponding excited half-integer-spin multiplets. The weak current coupling 
of these fictitious hadrons is fixed by the corresponding transition probabilities near zero 
recoil. Namely, the (unpolarized) zero-recoil structure functions are expressed through the 
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effective transition amplitudes; on the other hand, they are given by the 1/rriQ expansion, 
in our case to the second order, of the actual 5-meson zero-recoil semileptonic structure 
functions, cf. Eqs. (g, flA~28j ). 

There are more structure functions for a particle with spin than for actual B mesons. 
However, considering the unpolarized states (for instance, averaging over spin) they are 
reduced to the standard ones, and we use the same notations for them for decays of f2 
as for B mesons in QCD. V-A interference (which is not relevant here) would require a 
nonvanishing recoil kinematics. 

The unpolarized structure functions for the transitions into the ^ + states are the same 
as the tree-level ones for B decays [4*7] : 

t»i=9l (mi+ ? )2 ~ g2 +9 2 v {mi ~ m ) 2)2 ~ q \ w 2 = 2(g 2 A +g 2 v )ml w 3 = 2g A g v m 1 , (A.29) 



with the overall factor JE 25{qo — mi I"' 2+q ), where we have used the notation of Eq. (86), 



"> 2 — TO?- 1 -" 2 

2m 

namely 

(| + |cT7i75&|fio) = 5-aX7m75*o, (§ + |cT7i 6 l n o) = flVrV^o ■ (A.30) 

Note that here we use the full bispinors and assume their relativistic normalization. In this 
context the mass m 1 refers to the ground state in the beauty sector, m 1 ~m b +A while m 2 
to the excited state for charm, m2 — m c +A+e ra d- Only W\ } 2 contribute to the total width; 
it is given in Eqs. (87), (88). 

For transitions into | radial excitations we employ 



J + |c7 At 7 5 6|fi ) = 9aX^o, (§ + |c7^|fi ) = SVX/^75*o, (A.31) 



where Xn are likewise fully relativistic Rarita-Schwinger wavef unctions. (The vector matrix 
element is considered only for completeness; its transition amplitudes into | + states involve 
further suppression: either more powers of velocity, or extra 1 / vtlq or an additional overall 
a s .) With this convention we have 



(m 1 +m 2 ) 2 -q 2 a (m 1 -m 2 ) 2 -q 2 \ ( m 2 m 2 , 2 m \ 
3^ +9v 3^ j^--^), w 2 = ^ 



w 3 = 0, W4 = m 2 wi, w 5 = ~Wi. (A. 32) 

Integrating over the full phase space we obtain the corresponding width 



,(1) _ G\Ml\V ch \ 2 , (|) , , M l 



T £v- 792^3 9a,v z A,v( r )> r_ Mf' (A.33) 
with the phase space factors 
z%(r) = \(z { }\r)±~z^\r)\ zf' \r) = — - 2r - § r 2 + 3r 3 - §r 4 -4r 2 lnr, 



6a/t 



4-r 3 / 2 (12-f r-|r 2 + (6+8r)lnr) (A.34) 



66 



The singularity at r — >■ reflects here the ultraviolet problems of point-like higher-spin 
particles. Note that all the components of the vector transition amplitude vanish at zero 
recoil; the zero-recoil analysis does not constrain gy for | + . However, for this very reason 
it is generally suppressed by a higher power of the SV parameter A/M: 

Ml (§) A 5 Ml a), , 11 A 2 A 5 



z\*\r) ~ — — , z)?'{r) ~ —— at A < M . 



192tt 3 y ! ~ 30vr 3 ' 192tt 3 y v ; " 84 M 2 30tt 3 

Numerically, the yield in the vector channel is negligible for any relevant ratio M3/M0; the 

vector current does not contribute to the production of | + in our approximation. 

The effective formfactors at zero recoil are obtained comparing the structure functions 



in Eqs. (A. 29) or ( A. 32 ) in the zero-recoil kinematics (which is q 2 = q 2 , = (mi — rri2) 2 ) with 
those in Eqs. (52~j), (A. 28). In this way the constants g^, gy for | + and g& for | + above are 



expressed through the amplitudes introduced in Sec. 6.1 



£ 9a ] - + fi> £9v ] - (^-^) ^A-^/h^ 9 ' ( A - 35 ) 

where e is the mass gap for a particular state. 
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